AN  INVESTIGATION  OF  LOW  LAMBDA 
YIELD  FAILURE  DESIGN  OF  STIFFENED 
CYLINDRICAL  PRESSURE  VESSELS  UNDER 
EXTERNAL  HYDROSTATIC  PRESSURE 


Daniel  George  Hickey 


AN  INVESTIGATION  OF  LOW  LAMBDA 
YIELD  FAILURE  DESIGN  OF  STIFFENED 
CYLINDRICAL  PRESSURE  VESSELS  UNDER 

EXTERNAL  HYDROSTATIC  PRESSURE 

by 

DANIEL  GEORGE  HICKEY 

ENSIGN,  UNITED  STATES  NAVY 

B.S.,  United  States  Naval  Academy  (1971) 

SUBMITTED  IN  PARTIAL  FULFILLMENT 
OF  THE  REQUIREMENTS  FOR  THE 
DEGREE  OF  MASTER  OF  SCIENCE 
IN  NAVAL  ARCHITECTURE  AND  MARINE  ENGINEERING 

at  the 
MASSACHUSETTS  INSTITUTE  OF  TECHNOLOGY 
June,  1972 


/  if-t 


AN  INVESTIGATION  OF  LOW  LAMBDA 
YIELD  FAILURE  DESIGN  OF  STIFFENED 
CYCLINDRICAL  PRESSURE  VESSELS  UNDER 

EXTERNAL  HYDROSTATIC  PRESSURE 

by 
Daniel  G.  Hickey 

Submitted  to  the  Department  of  Naval  Architecture  and 
Marine  Engineering  in  partial  fulfillment  of  the 
requirements  for  the  degree  of  Master  of  Science  in 
Naval  Architecture  and  Marine  Engineering. 

ABSTRACT 


The  designer  of  stiffened  cylindrical  pressure  vessels 
is  often  faced  with  the  problem  of  comparing  structures 
of  wiedly  varying  geometries.   To  do  this,  he  makes  use 
of  a  number  of  dimensionless  parameters  whose  function  is 
to  reduce  the  number  of  variables  in  the  design  problem. 

The  object  of  this  thesis  was  to  investigate  the 
effect  of  the  dimensionless  parameter,  X,  on  stiffened 
cylindrical  pressure  vessel  design.   In  particular,  the  A 
region  where  collapse  of  vessels  due  to  yield-initiated 
buckling  between  the  frames  was  identified  and  delimited. 
This  region  was  used  as  the  first  iteration  in  a  design 
procedure  aimed  at  producing  a  vessel  of  maximum  structural 
efficiency.   Design  curves  which  may  be  used  to  develop 
such  a  vessel  are  provided  for  shell  thicknesses  up  to 
h/Dm  =  .01.   These  curves  are  based  on  both  maximum 
principle  stress  and  Hencky-von  Mises  stress  criteria  and 
provision  is  made  in  a  computer  program  for  use  of  other 
methods.   A  modified  Foppl  Formula  is  used  for  frame 
evaluations. 

The  formulas  used  range  from  conservative  to  liberal 
in  their  approximations  of  collapse  pressures.   However, 
use  of  even  the  most  conservative  formulas  results  in 
structural  efficiencies  which  compare  favorably  with 
present-day  practices. 
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NOMENCLATURE 

A    Cross  sectional  area  of  frame,  in 

Afc>    Cross  sectional  area  of  frame  plus  plating  of 
width  b,  in2 

Ai/A2,A3   Plasticity  coefficients,  see  Reference  (10) 

b    Width  of  frame  at  shell,  in 

c     Distance  from  neutral  axis  of  bending  to  extreme 
outside  fiber,  in         __ 

D    Outside  diameter  of  shell,  in 

d    Depth  of  rectangular  frame,  in 

D]-,    Diameter  to  neutral  axis  of  frame  plus  plating  of 
width  b,  in 

Dm  Diameter  to  neutral  axis  of  shell  alone,  in 

E  Modulus  of  elasticity,  lbs/in 

e  Eccentricity,  in 

Es  Secant  modulus  of  elasticity,  lbs/in 

f  Width  of  flange  of  T  frame,  in 

h  Shell  thickness,  in 

Ib   Moment  of  inertia,  frame  plus  plating  of  width  b, 
in1* 

Ie   Moment  of  inertia,  frame  plus  plating  of  effective 
width,  see  Reference  (10),  in 

If  Moment  of  inertia,  frame  alone,  in"* 

I^  Moment  of  inertia,  frame  plus  plating  of  width  L,  in 

L  Unsupported  span  of  plating,  in 

L^  Bulkhead  spacing,  in 

Lf  Frame  spacing,  in 

L-,  Unsupported  span  of  plating  taken  as  frame  sapcing,  in 


n    Number  of  complete  waves  in  buckled  configuration  - 
circumf erentially 

pc    Collapse  pressure,  lbs/in 

p-j_    Actual  collapse  pressure  desired,  lbs/in2 

Pcl'Pc2/Pc3   Temporary  values  of  pc  for  iteration 

qf    Collapse  load  per  unit  axial  length  of  a  ring  frame, 
lbs/in2 

R    Outside  shell  radius,  in 

R^    Radius  to  neutral  axis  of  "frame,  in 

Rm    Radius  to  neutral  axis  of  shell,  in 

R]_    Radius  to  neutral  axis  of  frame  plus  shell  of  width 
L,  in 

V    Volume  displaced  by  cylinder  per  unit  length,  in3/ft 

Wt   Weight  per  unit  length,  lbs/ft 

A    Weight  of  salt  water  displaced  by  cylinder  per  unit 
length,  lbs/ft 

n  Efficiency  -  defined  in  Section  III 

X  Slenderness  Ratio  -  defined  in  Section  I 

u  Poissons  ratio 

of  Axial  stress  at  mid-thickness,  lbs/in2 

Oip    Tangential  (circumferential)  stress  at  mid-thickness, 
lbs/in2 

Oy        Material  yield  stress,  lbs/in2 

ij>     Pressure  factor  -  defined  in  Section  I 
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I .   INTRODUCTION 

Stiffened  cylindrical  pressure  vessels  under  external 
hydrostatic  pressure  fail  by  one  or  any  combination  of 
three  collapse  modes.   While  both  theoretical  and 
empirical  expressions  have  been  developed  to  determine 
collapse  pressures  for  all  modes,  it  is  up  to  the  designer 
to  determine  what  structural  proportions  will  result  in 
the  most  efficient  use  of  material  strength  and  stiffness 
properties  for  a  required  collapse  pressure.   Two  dimen- 
sionless  parameters  which  have  been  widely  used  to  define 
different  geometries  and  their  capabilities  are   iJj  and 
X.  \\)    is  called  pressure  factor  and  it  is  a  ratio  of  a 
structure's  actual  collapse  pressure  to  the  collapse 
pressure  of  an  unstiffened  tube  with  the  same  h/Dm  ratio. 


H>      = 


m  J 


X  is  called  slenderness  ratio.   It  is  analogous  to  the 
slenderness  ratio,  1/r,  of  column  theory,  however,  it 
includes  terms  which  allow  for  different  material  properties. 


X   = 


'VD^ 


(h/Dm) 3 


°Y 


Despite'  t  .ieir  fairly  wide  use  by  designers  in 
comparisons  of  different  structural  geometries,  this 
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researcher  is  unaware  that  any  attempt  has  been  made 
to  develop  a  relationship  between  these  parameters  which 
would  indicate  their  effect  on  structural  efficiency. 
Also,  it  has  been  stated  that  these  parameters  are 
sometimes  misleading  (7) ,  however,  no  definite  limits 
on  their  usefulness  have  been  established  in  the  litera- 
ture, and  only  a  few  other  dimensionless  parameters  have 
been  developed  to  be  used  instead.   A  need  exists  therefore, 
to  define  the  limits  of  usefulness  of  both  \\>   and  X  and 
to  develop  within  those  limits,  relationships  between 
them  and  structural  efficiency  which  may  be  readily 
applied  in  the  early  design  stages  of  a  stiffened  cylindrical 
shell. 

In  this  report,  two  failure  criteria,  maximum  circum- 
ferential stress  and  Hencky-von  Mises,  have  been  used 
to  develop  relationships  between  X,  ty ,      and  efficiency. 
In  the  course  of  the  development,  certain  trends  of  the 
parameters  have  been  noted  and  in  the  concluding  remarks, 
useful  limits  of  X  and  ty   are  indicated.   In  particular, 
the  "low  lambda"  region  has  been  defined  and  an  intro- 
duction made  to  the  concept  of  maximizing  structural 
efficiency  within  that  region.   Design  curves  are  presented 
to  enable  optimization  within  the  low  lambda  region  of 
structures  wi+-h  collapse  pressures  up  to  2500  psi  and  examples 
of  the  process  are  given.   Special  note  should  be  taken 
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of  Appendices  B  and  C  where  remarks  concerning  very 
low  lambda  design  are  made. 
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II.   DEVELOPMENT  OF  DESIGN  CRITERIA 

Historically,  the  largest  and  most  critically 
designed  stiffened  cylindrical  pressure  vessels  have 
been  intended  for  use  as  submersibles .   Such  vessels, 
under  external  hydrostatic  pressure,  have  been  known  to 
collapse  in  one,  or  a  combination  of  three  modes.   Since 
the  mechanism  which  precipitates  failure  in  each  mode 
is  different,  the  designer  must  take  all  three  modes 
into  consideration  when  designing  for  a  given  collapse 
pressure.   However,  he  has  a  choice  in  the  design  method 
he  may  use  when  developing  a  highly  efficient  structure. 
He  may  design  toward  one  mode  of  collapse,  seeking  to  eli- 
minate the  possibility  of  collapse  in  others  while  changing 
structural  proportions  to  obtain  a  high  efficiency.   His 
alternative  is  to  simply  design  toward  the  most  efficient 
structure,  not  worrying  about  the  mode  of  collapse  as  long 
as  the  required  collapse  pressure  is  obtained. 

An  example  of  the  first  method  may  be  found  in  the 
design  of  the  relatively  shallow-diving  submarines  of 
World  V7ar  II.   These  vessels  were  designed  to  collapse  by 
asymmetric  buckling  between  frames.   It  was  felt  that 
imminent  collapse  would  be  signalled  by  the  appearance 
of  lobes  girdling  the  hull  and  the  crew  would  be  given  a 
chance  to  surface  before  total  collapse  occurred  (22) . 
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This  explanation  indicates  that  for  the  volume-limited 
submarines  of  that  era,  efficiency  criteria  could  be 
established  to  a  large  degree  by  safety  considerations. 

For  deeper  diving,  weight-limited  submarines, 
the  first  method  proved  inadequate  and  the  second  method 
was  refined  in  a  number  of  ways  to  allow  development  of 
highly  efficient  vessels.   A  refinement  which  will  now 
be  discussed  is  one  which  assumed  that  ty    is  limited  to 
1.0.   This  assumption  used  the  maximum  circumferential 
stress  criterion  and  minimized  weight  by  using  the 
smallest  possible  number  of  frames.   It  resulted  in  a 
structure  proportioned  such  that  X  =    .8.   In  order  to 
visualize  this  approach,  figure  1  has  been  developed. 

Figure  1  is  a  plot  of  ty   versus  d/h  (frame  depth/ 
shell  thickness)  for  four  X  values.   As  implied  by 
the  abscissa  scale  which  indicates  relative  frame  strength, 
the  frames  are  rectangular  and  only  their  depth  is  varied 
to  vary  strength.   The  basic  geometries  used  in  this  fig- 
ure are  taken  from  a  report  by  Trilling  (19)  and  one  of 
his  curves  (curve  2) ,  as  well  as  his  experimental  data 
(crosses)  are  plotted.   The  other  three  curves  plotted 
on  figure  1  represent  the  effect  of  changing  X.   In  all 
three  cases,-  J  he  distance  between  frames  was  the  only 
parameter  varied.   Values  of  ty   were  obtained  by  using 
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von  Sanden  and  Gunther ' s  formula  #92A  to  predict  collapse 
pressure. 

It  will  be  noted  that  curve  3  of  this  figure  lies 
almost  entirely  along  the  line,  \\>   =   1.0.   The  X  value 
of  this  curve  is  .814,  and  it  is  this  relationship  which 
forms  the  basis  of  the  design  method  cited  above.   Just 
as  the  earlier  method  of  designing  for  a  particular 
collapse  mode  proved  inadequate,  this  method  has  also 
been  unable  to  develop  the  efficient  structures  necessary 
to  resist  high  pressures.   The  remainder  of  this  section 
will  be  devoted  to  developing  a  design  method  which  will 
result  in  highly  efficient  structures.   To  begin  this 
development,  a  further  look  at  figure  1  is  necessary. 

On  the  basis  of  curve  2  and  the  experimental  points, 
one  of  Trilling's  conclusions  was  that  "except  for  very 
short  thick  tubes,  .  .  .  according  to  formula  #92A  the 
stronger  the  frame,  the  weaker  the  vessel."   This  was 
explained  by  noting  that  very  strong  frames  form  "hard 
spots"  when,  under  external  pressure,  they  contract  less 
than  the  shell.   The  hard  spots,  in  turn,  cause  stress 
concentrations  in  the  otherwise  uniformly  contracting 
shell  and  lower  collapse  pressures. 

Curve  4  indicates,  however,  that  this  conclusion  is 
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not  necessarily  true,  even  for  thin  shells  and  large  frames. 
The  only  difference  between  the  model  of  curve  4  and 
Trilling's  models  is  the  decreased  frame  spacing  which 
results  in  different  A  values;  .898  for  Tilling' s  models 
(curve  2)  and  .68  3  for  curve  4.   According  to  curve  4, 
shell  strength  actually  increases  with  frame  size  when  the 
frames  are  relatively  close  together.   The  seeming  contra- 
diction of  Trilling  may  be  explained  by  an  extension  of 
his  logic.   If  the  frames  causing  stress  concentrations 
are  closer  together  (low  lambda)  the  overall  stress 
pattern  will  be  uniform  enough  that  premature  yielding 
and  buckling  will  not  take  place. 

The  obvious  lower  limit  of  X  is  obtained  when  the 
frames  are  in  contact  (X-*0)  (24)  thereby  creating  a  thicker 
shell  (rectangular  frames)  or  a  webbed  sandwich-type 
construction  (T-frames) .   Both  of  these  geometries  result 
in  putting  the  frame-caused  stress  concentrations  in  such 
close  proximity  that  the  overall  stress  pattern  approaches 
uniformity.   This  uniformity  of  stress  pattern  implies  that 
the  yield  strength  of  most  of  the  material  in  the  structure 
will  be  utilized  before  collapse  occurs,  which,  in  turn, 
implies  that  a  structure  of  very  high  efficiency  is  possible 
in  this  X  area. 

With  a  lower  limit  of  X  established  at  0,  an  upper 
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limit  is  necessary  to  define  the  range  of  geometries 
where  the  effect  of  increasing  frame  strength  is  not 
detrimental  to  collapse  pressure.   This  limit  will  be 
referred  to  as  Xm.   For  structures  with  X  values  above 
Xm,  increasing  the  frame  strength  weakens  the  vessel. 
For  X  values  below  X  ,  increasing  the  frame  strength 
strenthens  the  vessel.   An  initial  value  of  this  limit 
can  be  found  by  direct  use  of  Windenburg ' s  approximation 
of  von  Mises1  formula  for  collapse  between  frames  due 
to  elastic  shell  instability.   The  development  of  this 
formula  is  outlined  in  Appendix  A.   In  its  simplest 
form,  this  formula  is: 

ip  =    1.3 


Recalling  that  tJj  =    pc     ,  Windenburg's  formula 

determines  instability  collapse  in  terms  of  the  hoop  stress 
relation,  i.e.  , 

.  ho,,   1.3 

pc  =  2irmy  — r~ 

therfore,  when   1.3    =  1.0,  failure  is  due  primarily 

X2 
to  yielding  of  the  material  rather  than  plastic  instability, 

At  this  point,  X  =  1.14.   If  X   1.14,  then  instability 

collapse  is  indicated  since  \\)      1.0.   In  this  situation, 
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strengthening  the  frames  cannot  delay  collapse.   It  may 

be  tentatively  assumed  that  when  X   1.14,  strengthening 

the  frames  increases  overall  strength.   If  this  is  so, 

then  an  initial  approximation  of  X   is  1.14.   It  is  important 

to  note  that  this  is  a  very  approximate  method  for 

determining  Xm.   A  glance  of  Trilling's  results  on  figure  1 

will  indicate  that  collapse  at  \p   values  below  1.0  has  occurred 

in  structures  with  X  values  as  "low  as  .89.   Also,  with 

this  method,  Xm  is  independent  of  any  individual  structure 

parameters.   Since  this  is  questionable,  another  means 

of  determining  Xm  is  necessary. 

Returning  to  curve  3  on  figure  1,  it  will  be  recalled 
that  for  the  X  value  used  in  that  curve,  4>   =    1 . 0  for  all 
frame  sizes.   As  long  as  the  frames  are  of  adequate 
strength  or  stiffness  to  carry  a  given  percentage  of  the 
shell  load,  their  size  has  no  bearing  on  the  strength  of 
the  cylinder.   Figure  2  elaborates  on  this  point.   The 
curves  of  figure  2  are  collapse  pressures,  determined 
by  #92A,  plotted  as  fy   against  values  of  X.   Each  curve 
represents  one  size  of  frame  and  all  structures  have  the 
same  h/D  ratio,  so  that  only  the  distance  between  frames 
has  been  varied  to  obtain  changing  X  values.   The  point 
at  which  all  curves  cross  is  at  ty =   1.0  and  represents 
the  afore-mentioned  Xm.   Unlike  the  approximate  Windenburg 

determination  of  this  point,  in  this  determination,  Xm 
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depends  on  the  value  of  h/D  as  well  as  the  fact  that 
ty  =  1.   This  dependence  is  borne  out  in  the  derivation 
of  an  expression  for  Xm  from  formula  #92A.   This  deri- 
vation is  found  in  Appendix  E  and  results  in  the 
equation: 


*m 


2.055Dm      uy  (A) 


h         E 

which  was  used  to  plot  h/Dm  versus  Xm  for  different 
values  of  ay  on  figure  3. 

Figure  3  may  be  used  to  determine  for  a  given 
h/D   and  a   what  X  is  necessary  to  prevent  collapse  be- 
tween the  frames  due  to  instability  (i>   <r  1.0)  .   Areas 
to  the  right  of  the  curve  in  question  are  areas  of 
instability  collapse,  those  to  the  left,  primarily 
yield  failure.   As  an  example  of  this,  Trilling's 
experimental  results  have  again  been  plotted.   Due  to 
slightly  varying  geometries  and  a    ' s ,  all  the  models  he 
used  fall  within  the  areas  of  the  dots  shown.   The 
figure  indicates  that  for  a   =  30,000  psi  (see  appendix) 
all  the  models  failed  at  a  \\>   below  1.0.   Reference  to 
figure  1  shows  that  this  was  indeed  the  case, and  a 
further  glance  at  Reference  (19)  indicates  that  a  number 
of  Trilling's  models  failed  by  lobar  buckling  or  general 
instability,  rather  than  by  yield  initiated  collapse. 
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In  Appendix  A,  a  number  of  formulas  and  yield 
criterion  are  cited  as  being  more  accurate  than  the 
maximum  circumferential  stress  criterion  in  determining 
collapse  pressure  due  to  yielding.   One  of  these,  Wenk '  s 
formula,  has  been  investigated  using  the  same  models 
used  in  figures  1-3  in  order  to  ascertain  the  influence 
of  increased  accuracy  on  the  computation  of  Xm.   Figure 
4  is  the  result  of  this  investigation. 

Wenk  applied  formulas  #92  and  #92A  to  the  Hencky- 
von  Mises  yield  criterion  at  mid-thickness,  mid-bay, 
and  obtained  collapse  pressures  within  nine  percent  of 
actual  experimental  values.   The  Hencky-von  Mises  criterion 
in  two  dimensions  is: 

Oy  =   arp2  -  OrpOp   +   aF2 

When  it  is  applied  to  an  unstiffened  shell  where  aT  = 
2aF,  yield  failure  occurs  at  aT  =  1.16av,  i.e.,  x/j   = 
1.16.   Just  as  figure  2  indicated  that  for  the  maximum 
circumferential  stress  criterion,  there  is  a  X  value 
where  frame  size  has  no  effect  on  cylinder  strength,  figure 
4  shows  that  the  same  is  true  for  the  Hencky-von  Mises 
criterion.   However,  the  X  value  in  this  case  occurs  when 
ty   =    1.16  instead  of  1.0.   Wenk's  formula: 
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F2  -  .5F  +  .25 


like  formula  #92A,  is  actually  the  hoop  stress  relation 
times  a  factor (defined  in  Appendix  E) .  When  it  is  solved 
for  F  by  setting  \p   =   1.16,  a  Xm  versus  h/D   relation 
similar  to  formula  (A)  may  be  found.   This  formula  is 
shown  below  and  a  detailed  derivation  of  it  is  given  in 
Appendix  E: 


2-65  Dm       *y  (B) 


Xm 


If  plotted  in  the  form  of  figure  3,  this  equation  would 
result  in  a  shift  of  all  curves  . 11A  to  the  right  of 
those  obtained  from  #92A.   This  shift  reflects  the  conser- 
vative nature  of  von  Sanden  and  Gunther's  formula 
and  criterion. 

Thus  far,  equations  and  criteria  have  been  developed 
to  determine  the  extreme  values  of  X   between  which 
stiffened  cylindrical  shell  failure  is  primarily  due  to 
yielding  of  the  material  between  the  frames.   Justification 
for  using  this  as  the  design  collapse  mode  comes  from  a 
number  of  sources.   Wenk,  a  major  advocate  of  the  "one- 
hoss-shay"  collapse  philosophy,  states  that  the  approach 
of  designing  for  yield  failure  is  not  inconsistent  with 
that  philosophy  (221.  Eecause  of  the  effects  of  such 
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unknowns  as  material  residual  stresses  and  out  of  roundness, 
design  for  collapse  by  all  three  modes  simultaneously 
is  difficult,  if  not  practically  impossible.   Therefore, 
the  mode  which  is  least  affected  by  such  strength  modifiers 
should  be  the  design  collapse  mode.   Conversely,  the 
collapse  mode  most  affected  by  the  modifiers  should  be  the 
first  mode  taken  into  consideration  and  designed  out  of 
existance.   Out  of  roundness  is-  the  primary  reason  for 
local  failure  of  frames  which  can  precipitate  premature 
general  instability  collapse  (15)  .   For  this  reason, 
collapse  by  general  instability  has  been  noted  as  being 
of  prime  concern  in  the  preliminary  design  (10) . 

Reference  to  formulas  for  general  instability  given 
in  the  appendix  will  reveal  that  this  collapse  mode,  unlike 
the  others,  is  a  function  of  cylinder  length  (between 
bulkheads) ,  as  well  as  frame  area,  moment  of  inertia,  and 
spacing.   Judicious  selection  of  bulkhead  spacing  (if  a 
free  variable)  and  minimum  frame  sizes  can  eliminate 
this  collapse  mode  from  consideration.   This  should  be 
done  in  the  first  stages  of  design  development.   In  the 
final  stages  collapse  pressure  for  this  mode  should  again 
be  checked  and  out  of  roundness  taken  into  consideration. 

V7ith  general  instability  designed  for,  collapse 
between  the  frames  due  to  buckling  should  be  eliminated 
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by  insuring  that  the  final  design  X  will  be  between  0 
and  Xm.   This  means  that  the  yield  strength  of  the  material 
will  be  utilized  as  fully  as  possible  by  insuring  against 
premature  failure  initiated  by  instability.   In  this  way, 
the  benefits  of  using  high  strength  materials  will  not 
be  nullified  by  instability  collapse  at  low  pressures. 

In  the  next  section,  considerations  for  the  final 
step  of  designing  for  yield  initiated  failure  between  the 
frames  will  be  developed.   First,  however,  the  question 
of  frame  adequacy  must  be  looked  into. 

Figures  1,  2,  3,  and  4  were  all  developed  from 
analysis  of  cylindrical  shells  with  rectangular  stiffeners. 
T-shaped  stiffeners  are  in  general  use  today  however,  so 
the  question  may  be  raised  as  to  what  effect  the  shape 
of  the  stiffeners  has.   An  investigation  of  the  formulas 
for  yield  failure  previously  cited  will  reveal  that  they 
are  all  based  on  the  assumption  of  adequate  frame  strength. 
The  only  frame  parameters  employed  in  them  are  flange 
width  in  contact  with  the  shell  and  frame  area.   At   X  - 
Xm,  even  these  parameters  cease  to  influence  collapse 
pressure.   Therefore,  it  can  be  stated  that  frame  shape 
has  no  influence  on  formulas  (A)  or  (B) .   In  general, 
as  long  as  adequate  strength  is  maintained,  frame  shape 
does  not  influence  the  collapse  pressures  of  stiffened 
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cylindrical  shells. 

The  question  of  what  constitutes  adequate  frame 
strength  has  been  directly  approached  by  a  number  of 
investigators  (2,11,19).   The  standard  practice  of 
von  Sanden  and  Gunther  (20)  and  other  early  designers 
was  to  design  the  frame  to  be  able  to  withstand  1.1 
times  the  total  pressure  on  the  length  of  plating  it 
supported.   This  procedure,  which  employed  the  Foppl 
formula,  was  found  to  result  in  "excessively  heavy  frames" 
(19) .   In  addition,  the  Foppl  formula  has  been  determined 
to  give  conservative  results  for  frame  buckling  pressures. 
In  light  of  these  facts,  the  frame  design  requirement 
chosen  for  the  purposes  of  this  report  was  that  the  frame 
be  strong  enough  to  support  no  more  than  the  load  on 
one  effective  length  of  shell.   The  conservativeness  of 
the  Foppl  formula  was  accounted  for  by  assuming  that  it 
predicted  frames  alternatively  4,  5,  or  6  times  too 
strong.   The  resulting  formula  used  is: 

24EIb 

P„  frame    =   p_  shell   =   X  (4, 5, or  6) 

D>b  Li 

Before  proceeding  into  the  actual  design  iteration 
resulting  from  the  establishment  of  design  criteria  for 
stiffened  cylindrical  pressure  vessels  under  hydrostatic 
pressure,  a  recapitulation  of  the  criteria  is  in  order: 
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1.  Upon  definition  of  -operating  requirements,  initial 
Structure  parameters  should  be  established  to 
minimize  the  possibility  of  general  instability 
collapse  of  frame  and  shell. 

2.  Possibility  of  collapse  due  to  instability 
between  frames  should  be  reduced  by  limiting 
design  choices  to  geometries  located  between 
A  =  0  and  X   =   \m. 
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III.   VESSEL  EFFICIENCY 

The  measure  of  efficiency  used  in  this  paper 

was  suggested  by  Wenk  in  (22) .   It  is  slightly  more 

complicated  than  the  traditional  Weight/Displacement 

ratio  in  that  it  takes  into  account  collapse  pressure. 

p  V 
^c 


105W 

pc   is  the  collapse  pressure  -  lbs/in2 

V    is  the  volume  per  unit  length  of  cylinder  -  in3/ft  length 

W    is  weight  per  unit  length  of  cylinder      -  lbs/ft  length 

When  used  in  a  qualitative  manner,  this  term  clearly 
shows  the  relative  efficiencies  of  different  geometries. 
However,  it  is  not  independent  of  material  properties 
such  as  density  and  yield  stress.   The  efficiencies  in 
this  report  have  been  calculated  using  HY  80  as  the 
material.   In  order  to  apply  them  to  materials  of  other 
densities  and  yield  strengths,  they  must  be  multiplied 
by  the  following  factors: 

oy  ,,  .284 

solooo   and/or    —p- 

For  example,  a  geometry  which  would  have  an  efficiency 

of  1.25  for  HY  80,  has  an  efficiency  of: 

1  25a50,000  ,.284    =   2  ig 
l'Zb    80,000   ~        ^•iy 
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when  the  material  is  50,000  psi  aluminum.   Once  found  in 
the  above  manner,  the  efficiency  may  be  used  to  calculate 
the  traditional  Wt/Displacement  ratio  using  the  following 
conversion  formula: 

Wt    _   J^c 
A       3720n 

Reasons  for  this  manipulation  of  n  will  become  more  clear 
at  the  end  of  the  next  section. 
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IV.   DESIGN  ITERATION 

In  the  preceding  section,  design  criteria  have 
been  developed  which,  if  used,  will  produce  a  generalized 
cylinder  geometry  whose  A  value  will  be  between  0  and  Xm. 
This  will  result  in  a  collapse  pressure  factor  in  excess 
of  1.0,  thus  insuring  that  collapse  will  be  primarily 
due  to  yield  and  the  yield  strength  of  the  material  will 
be  highly  utilized.   This  section  will  outline  development 
of  an  iterative  method  to  find  the  particular  X   and  i|> 
values  which  result  in  maximum  structural  efficiency 
for  a  given  collapse  pressure. 

Figures  5,  6,  and  7  are  plots  of  efficiency  against 
X  for  different  values  of  h/Dm  and  frame  strength  factors. 
As  stated  previously,  Foppl's  frame  formula  has  been  used 
in  a  modified  form  which  arbitrarily  reduces  its  apparent 
conservativeness .   In  determining  frame  size  for  figure 

5,  Foppl's  formula  was  multiplied  by  a  factor  of  4 .   In 
figure  6,  the  frame  factor  was  5  and  in  figure  7,  it  was 

6.  The  effect  of  this  increasing  theoretical  strength 
of  frames  was  to  reduce  frame  size  and  increase  apparent 
structure  efficiency  in  succeeding  figures.   Since  this 
is  an  arbitrary  change,  model  tests  are  necessary  to 
determine  the  actual  frame  factor  which  should  be  used. 
It  was  felt  that  a  range  of  factors  such  as  that  covered 
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would  indicate  design  trends  which  may  be  beneficial  in 
general  use  without  the  necessity  of  establishing  rigorous 
theoretical  frame  design  criteria.   Further  substantiation 
of  this  method  will  be  provided  later  in  this  section. 

The  formula  used  to  determine  collapse  pressures 
for  these  three  figures  was  #92A  (maximum  circumferential* 
stress  criterion) .   Although  this  is  now  considered  to 
be  less  accurate  than  a  number  of  possible  formulas,  it 
has  the  advantages  of  always  being  conservative  and 
having  been  widely  used.   The  curves  in  figures  5,  6,  and 
7  indicate  that  for  a  given  h/Dm  and  frame  factor,  the 
efficiency  increases  with  X  until  a  maximum,  nmax/  is 
reached.   Thereafter,  the  efficiency  steadily  decreases, 
with  increasing  X  .   It  is  this  point,  nmax/  that  the 
designer  of  stiffened  cylindrical  pressure  vessels  is 
interested  in,  therefore,  only  enough  points  have  been 
plotted  to  indicate  the  hump  of  each  efficiency  curve. 

An  interesting  trend  which  is  readily  apparent  in 
all  three  figures  is  that  the  efficiencies  of  all  h/D 
combinations  lie  on  the  same  line  at  the  point  where 
X  =  Xm   for  each  h/Dm  (see  figure  3).   In  figure  5,  this 
line  has  a  slope  of  about  -.5.   In  subsequent  figures,  the 
slopes  appear  the  same,  so,  since  the  lines  begin  at  higher 
values,  the  final  efficiencies  will  be  higher.   Of  course, 
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extension  of  the  lines  upward  to  maximum  efficiency  also 
requires  extension  to  the  left  to  X  =  0  and  extension  of 
the  series  of  h/D  curves  to  extremely  thick  cylinders. 
This  extrapolation  of  computed  data  can  be  justified  on 
the  grounds  that  only  very  short,  thick  cylinders  fail 
by  pure  yielding.   Failure  in  this  way  takes  the  fullest 
possible  advantage  of  material  yield  strength  and  results 
in  the  highest  possible  structural  efficiency.   Since 
efficiency  in  this  report  is  also  a  function  of  displace- 

» 

ment  to  weight  ratio,  the  veracity  of  that  previous 
statement,  in  the  general  sense,  is  dependent  on  material 
density. 

Another  point  which  is  illustrated  in  figures  5-7, 
is  that  for  each  h/Dm  there  is  only  one  value  of  A  at 
which  maximum  efficincy  is  reached,  regardless  of  frame 
factor.   The  effect  of  the  frame  factor  is  to  change  the 
actual  maximum  efficiency  reached,  but  it  has  no  effect 
on  the  X  value  where  it  is  reached.   This  direct  relation- 
ship between  h/D   and  X  at  maximum  efficiency  is  shown 
in  figure  8.   It  may  be  explained  in  part  by  the  fact  that 
neither  X  nor  h/Dm  contain  any  frame  size  parameters  and 
both  Wenk's  formula  and  formula  #92A  assume  adequate 
frame  strength.   Also  at  that  point  of  maximum  efficiency, 
there  are  dis.  rete  values  of  \p   which  are  dependent  on 
frame  factor.   In  order  to  show  the  relationship  of  frame 
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f actor,  A  and  i»   at  maximum  efficiency,  figures  9  and 
10  have  been  plotted. 

Before  entering  into  the  use  of  these  figures, 
however,  a  slight  digression  is  necessary.   Recalling 
that  formula  #92A,  used  to  obtain  figures  5-7,  was 
shown  to  yield  less  accurate  collapse  pressures  than 
other  formulas,  it  became  necessary  to  investigate  a 
formula  of  higher  accuracy  to  determine  if  the  same 
trends  still  held.   The  formula  used  was  developed  by  Wenk 
and  has  been  introduced  in  Section  II.   Data  similar 
to  that  shown  in  figures  5,  6,  and  7,  was  developed  for 
the  same  geometries  using  Wenk ' s  formula.   The  actual 
results  are  given  in  Appendix  G.   Integrated  curves  from 
these  results  have  been  drawn  on  figures  8  and  10. 

Figure  8  is  a  graph  of  h/D  versus  A  for  points  of 
maximum  efficiency.   The  two  curves  represent  results 
of  the  use  of  formula  #92A  and  Wenk ' s  formula.   It  will 
be  noted  that  Wenk ' s  formula  shifts  the  curve  to  the  right. 
Since  this  expression  has  been  known  to  predict  collapse 
pressures  in  excess  of  actual  values,  it  may  be  expected 
that  the  curve  for  actual  structures  will  lie  between  the 
curves  shown.   No  mention  of  frame  size  is  made  on  this 
graph,  as  in  neither  formula  did  frame  strength  alter 
the  relationship  of  h/D  and  A   at  maximum  efficiency. 
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Figures  9  and  10  are  graphs  of  n    ,  X   and  ty   for 
#92A  and  Wenk's  formula  respectively.   The  three  lines 
in  each  case  indicate  the  three  frame  factors  used. 
Comparison  of  these  figures  shows  that  increased  accuracy 
does  not  change  the  trends  previously  noted.   In  fact, 
the  n  versus  A  curves  are  exactly  alike.   The  25  percent 
increase  in  efficiencies  in  figure  10  over  ffigure  9  is 
another  reflection  of  the  conservativeness  of  formula 
#92A. 

With  development  of  figures  8-10  complete,  an 
iterative  method  to  determine  geometries  of  maximum 
efficiency  using  these  figures  can  be  outlined.   Geometric 
parameters  which  can  be  found  are  h/Dm  and  X .   Given 
material  properties,  distance  between  stiff eners  may 
be  computed.   Actual  stiff ener  parameters  are  left  up  to 
the  discretion  of  the  designer  as  the  only  requirement  of 
the  formulas  used  in  this  procedure  is  that  the  stiffeners 
be  of  adequate  strength  and  stiffness.   For  comparison 
purposes,  the  frame  parameters  used  in  this  report  and 
those  frame  dimensional  ratios  coinciding  with  maximum 
efficiency  structures  are  given  in  Appendix  C.   To  begin 
the  iterative  procedure,  the  information  which  must  be 
specified  is: 

1.   type  of  analysis  to  be  used 

a.   maximum  principle  stress,  #92A 
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b.   Hencky-von  Mises  stress  at  mid-thickness, 
mid -bay,  #92  and  #92A 

2.  frame  factor  to  be  used:  4,  5,  or  6 

3.  ultimate  collapse  pressure  desired 

4.  material  properties,  ay  and  E 

(for  these  curves  E  =  3  X  10 7) 

5.  outside  diameter 

An  example  of  the  procedure  is  given  on  the  next  page, 
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Example  1: 

Given:        Use  of  #92A,  Frame  Factor  -  5 

cy   =   80,000  psi 
D    =   40  ft. 
p    =   1000  psi 


Steps: 


1.   Initial  h/Dm  guess  from  hoop  stress 
pc     =   2(h/D)oy      h/D  =  h/Dm 

(h/Dm)  =   1000/(2*80,000) 

(h/Dm)  =   .00692 

A      =   .48  from  figure  10 

i>  =1.29  from  figure  11 


=  1.29 


*v  °y 


.*.  pc  =  1290  psi  -  which  is  too  high 

Reduce  h/Dm  by  division 

Pel    =   100°    =  775  Psi 
1.29 

(h/Dm)  =   775/(2*80,000) 

(h/D  )  =   .00484 

'  m 

X      =   .58 

i>  =  i.4i 

p_     =   1.41*775  =  1090  psi  -  which  is 
still  too  high 
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3.   Reduce  h/Dm  again 

775        -in 

Pc2    =  T709~  =  71°  psl 

(h/Dm)  =   710/(2*80,000)  =  .00443 

X       =   .606 

i|i      =   1.44 

pc     =   1.44*710  =  1020  psi  -  which 

is  close  enough,  but  try  a  4th 

iteration 


4 .   Reduce  h/D 


m 


710        rQ_ 
Pc3     =  TT02        =   697  PS1 


(h/Dm)  =  697/(2*80,000)  =  .00436 

X  =  .613 

iJj  =  1.45 

pc  =  1.45*697  =  1012  psi 


Answer: 

X 

= 

.613 

h/D 

= 

.00436 

n 

= 

1.22 

* 

= 

1.45 

Results : 

D 

= 

40  ft. 

h 

= 

2.1  in. 

Lf 

- 

1.56  ft. 

Wt. 

= 

18,000  lb/ft 

conventional  V7t/Disp   =   .222 
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V.   DISCUSSION  OF  RESULTS 

The  results  of  this  report  may  be  divided  into  two 
classes;  those  due  to  analysis  and  those  due  to  synthesis. 
Analysis  by  formula  #92A  and  Wenk ' s  Formula  at  low 
lambda  values  has  shown  that  considerable  gains  in 
efficiency  can  be  realized  by  placing  frames  close  together 
and  reducing  shell  thickness.   This  results  in  i>   values 
greater  than  1  and  helps  to  realize  the  stress  potentials 
of  high  strength  materials.   The  extrapolation  of  the 
design  curves  in  figures  11  and  12  appears  to  justify 
the  use  of  sandwich  construction  for  very  high  pressures. 
Such  construction  is  analogous  to  design  in  the  very  low 
lambda  region,  and ' a  discussion  of  this  is  included  in  the 
appendix.   Inferences  which  may  be  drawn  from  the  treat- 
ment of  very  low  lambda  design  in  Appendix  B  concern  the 
apparent  inapplicability  of  X  and  ty  in  this  region.   This 
results  from  the  fact  that  in  this  region,  frame  strength 
contributes  so  much  to  overall  strength  that  it  cannot 
be  ignored  as  it  is  in  both  X   and  ty    . 

At  the  high  end  of  the  low  lambda  region,  maximum 
structural  efficiency  can  be  realized  only  in  relatively 
low  pressure  vessels.   This  is  shown  in  two  examples  in 
the  appendix. 
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The  result  of  synthesis  in  this  report  is  the 
iterative  method  outlined  in  Section  IV.   By  its  very 
nature,  this  method  results  in  thin-skinned  vessels 
with  frames  relatively  close  together.   The  questions 
which  this  geometry  may  raise  concerning  modes  of  collapse 
other  than  yield-initiated,  must  be  answered  by  indepen- 
dent analyses.   It  will  be  noted,  however,  that  the 
geometries  used  to  plot  the  curves  of  figures  8-13  were 
analyzed  for  failure  both  by  instability  between  frames 
and  general  instability  and  the  indicated  collapse  mode 
for  all  was  yield-initiated  failure  between  rings. 
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VI.   CONCLUSIONS 

1.  The  dimensionless  parameter,  A,  may  be  used  to 
define  all  stiffened  cylindrical  pressure  vessel  geometries 
except  those  in  which  the  frames  provide  little  support 
because  of  their  wide  spacing,  and  those  vessels  in  which 
the  frames  provide  most  of  the  strength  by  virtue  of  their 
size  and  close  spacing. 

2.  Absolute  structure  efficiencies  in  the  low  lambda 
region  are  highly  dependent  on  frame  spacing  and  only 
slightly  dependent  on  frame  size  (as  long  as  the  frames 
are  of  adequate  strength) . 

3.  Relative  structure  efficiencies  in  the  low  lambda 
region  are  virtually  independent  of  frame  size.  The  lambda 
values  at  which  maximum  efficiency  is  reached  depend  on 

the  h/Dn  ratio  and  are  not  affected  by  changes  in  frame  size, 

4.  Structures  with  \p   values  well  in  excess  of  1.0 
are  reasonable  and  generally  provide  increased  efficiencies 
over  those  where  \p    is  limited  to  1.0. 

5.  As  the  design  collapse  pressure  increases,  the 
value  of  lambda  which  coincides  with  maximum  structural 
efficiency  decreases.   This  implies  that  higher  efficiencies 
can  be  obtained  by  reducing  shell  thickness  well  below  that 
required  by  the  hoop  stress  relation  and  simultaneously 
increasing  the  number  of  frames. 

6.  At  extremely  low  lambda  values,  general  instability 
becomes  the  most  likely  mode  of  failure. 
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VII.   RECOMMENDATIONS 

Since  this  thesis  was  limited  to  studying  the 
parameter,  X,  and  its  effect  on  efficiency,  other 
definitive  parameters,  such  as  those  now  used  at  the 
David  Taylor  Model  Basin  should  be  investigated  for 
their  effect  on  efficiency.   It  has  been  noted  that  the 
lambda  term  is  relatively  useless  when  the  frames  are 
very  close  together  and  are  very  strong.   This  condition 
should  be  corrected  by  defining  a  new  parameter  similar 
to  lambda  in  its  application,  but  allowing  for  the  effect 
of  frame  strength  in  geometric  comparisons. 

The  basic  approach  of  this  thesis  should  be  used 
to  extend  the  curves  of  maximum  efficiencies,  found  in 
Section  IV,  into  the  regions  of  higher  collpase  pressures. 
With  this  extention,  a  thorough  re-examination  of  general 
instability  collapse  in  these  regions  should  be  undertaken, 
Finally,  experimental  tests  of  maximum  efficiency  geo- 
metries, as  developed  by  the  design  iteration,  should  be 
undertaken  to  provide  evidence  of  the  method's  worth. 
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APPENDIX  A 
HISTORICAL  BACKGROUND 

Prior  to  1917,  the  design  of  stiffened  cylindrical 
pressure  vessels  was  based  on  the  assumption  that  vessel 
strength  depended  primarily  on  frame  strength.   The  shell 
was  treated  as  a  flange  of  the  frame  and  strength  was 
calculated  from  equations  which  defined  collapse  by 
instability  of  frame  rings  under  radial  loads.   A  widely 
used  formula  for  this  was  developed  by  Levy,  but  is 
commonly  called  Foppl's  Formula  (19,24): 


3EI 

3f   = 


V 

While  this  formula  was  relatively  accurate  for  the 
time,  and  where  frames  alone  were  involved,  when  it  was 
applied  to  a  stiffened  cylindrical  shell,  its  shortcomings 
became  apparent  as  higher  collapse  pressures  were  required, 

In  order  to  improve  upon  this  frame-oriented  theory, 
von  Sanden  and  Gunther  derived  expressions  based  on  the 
observation  that  stiffened  cylindrical  shells  were  not 
collections  of  connected  rings,  but,  in  fact,  were  closed 
cylinders  reinforced  by  ring  frames.   Deciding  that  the 
shell  was  the  real  strength  member,  they  predicted  failure 
due  to  shell  yielding  at  the  area  of  maximum  stress.   This 
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maximum  principle  stress  criterion  could  be  satisfied 
by  the  longitudinal  stress  at  the  frames  or  by  the 
circumferential  stress  midway  between  them.   Formulas 
for  both  stresses  were  developed  and  the  well-known 
formulas  #92  and  #92A  resulted  when  pressures  were 
for  in  terms  of  the  stresses.   Substituting  material 
yield  stresses  into  #92  and  #92A  resulted  in  two  yield 
pressures  for  the  structure,  of  which  the  lower  one  was 
said  to  predict  failure.   It  must  be  emphasized  that 
failure  was  predicted  on  the  basis  of  the  onset  of 
yielding. 

For  a  number  of  years  after  the  1918  publication  of 
these  formulas,  #92  was  thought  to  give  the  critical 
collapse  pressure  since  the  greater  stress  was  longitudinal, 
In  1935,  however,  this  assumption  was  found  to  be  in 
error  when  circumferential  stress  was  determined  to  be 
critical  in  the  collapse  of  thin-shelled  pressure  vessels 
(19). 

The  von  Sanden  and  Gunther  formulas  are: 


#92  Py    =  2h/Dmay 


h+    1.8 IK  /.8  5  -  B^ 
I  1  +  6  ) 
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#92A     pv   =      2h/Dmgy 


y 


+  H  /.  85  -  B\ 
\    1  +   6J 


The  variables  K,  H,  B,  and  B  are  defined  in  the  Formula 
section. 

About  this  same  time,  it  was  generally  recognized 
that  the  area  between  stiffeners  very  rarely  collapsed 
as  a  result  of  pure  yielding.   In  anything  but  thick 
cylinders,  collapse  occurred  from  a  combination  of 
yielding  and  buckling  when  the  distance  between  stiffeners 
(effective  length)  was  below  a  critical  length.   This 
critical  length  was  dependent  on  shell  L/D  and  h/D  ratios 
as  well  as  material  properties.   When  the  distance  between 
stiffeners  was  above  the  critical  length,  the  tube  was 
considered  to  be  of  infinite  length  and  the  formula  of 
Breese  and  Bryan  for  collapse  by  instability  could  be 
applied: 


2E 

Pc    = ;    ft/Dm) 

1  -  V 


For  cylinders  below  the  critical  length,  formulas 
for  collapse  by  instability  were  developed  by  von  Mises 
and  Southwell.   Windenburg  simplified  von  Mises  lengthy 
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formula  and,  by  sacrificing . a  very  small  degree  of  accuracy, 
came  up  with  the  followina  formula  (25) : 


2.42E       (h/D)52 


(1  -  y2)3"   [L/Dra  -  .45(h/Dm)  ]12 

At  the  time  Windenburg  developed  the  above  formula, 
a  pc  versus  L/D  coordinate  system  was  commonly  used  to 
plot  experimental  and  theoretical  results  for  changing 
values  of  h/D.   In  an  effort  to  develop  an  analog  to  the 
slenderness  ratio,  L/r,  in  column  theory,  and  aid  comparison 
of  results,  he  developed  a  new  coordinate  system  based 
on  the  coordinates: 


Y  o 


Xc 

0 


2h/D 


m 


a/Dm)  ' 


(100h/Dm) 


Windenburg  then  corrected  his  coordinates  for  the  physical 
properties  of  materials  in  the  following  manner  and  came 
up  with  the  well-known  variabels,  slenderness  ratio  X  , 
and  pressure  factor,  ij>  . 

»o             Pc 
$  =        =  

oy  2h/Dmay 
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X  =     X 
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This  done,  Windenburg  rewrote  his  instability  formula  in 
terms  of  his  new  variables: 


* 


1.30 


X2-e 


3/t 

,  nAt.   1000a„    /  1   -   y2\  1 

where    e   =    .045  Y    I  \ 


.91 


lOOt/D 


may  be  neglected  in  general  to  get 


* 


1.30 


This  formula  is  plotted  on  figure  A-l  and  may  be 
compared  there  with  the  plot  of  von  Mises '  formula  and 
numerous  experimental  results.  One  thing  which  may  be 
seen  at  this  time  is  the  good  agreement  these  formulas 
have  with  eachother  and  the  experimental  results  above 
a  X  value  of  1.4. 


After  a  decade  of  relative  inactivity,  the  investi- 
gation of  stiffened  cylindrical  pressure  vessels  was  renewed 
in  the  early  1950' s.   The  main  intent  of  the  investigations 
of  this  period  was  to  correct  the  number  of  theoretical 
deficiencies  in  formulas  -92  and  #92A.   The  first  two 
deficiencies  were  that  the  formulas  predicted  linear 
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relationships  between  pressures  and  stresses  and  ignored 
the  effect  of  moments  developed  due  to  pressure  at  the 
ends  of  closed  cylinders.   Both  problems  were  corrected 
by  Salerno  and  Pulos  (17) ,  who  developed  stress  formulas 
similar  to  #92  and  #92A  but  taking  both  the  above 
moments  and  the  influence  of  axial  stress  in  the  shell 
tending  to  expand  the  frame  into  account.   Their  analysis 
resulted  in  a  non-linear  relationship  between  pressures 
and  stresses  which,  however,  did  not  alter  the  results 
of  stress  calculations  significantly  from  those  found  by 
#92  and  #92A. 

The  next  problem  found  with  von  Sanden  and  Gunther's 
theories  had  to  do  with  the  inadequacy  of  the  maximum 
principle  stress  criterion.   Use  of  the  Hencky-von  Mises. 
yield  criterion: 

a    2  *  aT2  -  a-pap  +  cf2 

at  the  exterior  mid-bay  point  resulted  in  predictions  closer 
to  actual  collapse  pressure  than  obtained  with  the  former 
criterion- (9).   However,  here  again  collapse  pressures 
were  associated  with  the  onset  of  yielding  and  no  allowance 
for  plastic  reserve  strength  was  made. 

This  condition  was  corrected  by  Wenk  et.  al.  (21)  in 
a  somewhat  haphazard  manner  by  taking  the  Hencky-von  Mises 
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yield  criterion  at  mid-bay  and  mid-thickness.   Although 
Vlenk   termed  this  a  "necessary,  if  not  sufficient  empirical 
yield  condition  to  produce  failure",  experimental  results 
have  since  indicated  that  thought  this  formula  predicts 
within  nine  percent  of  actual  collapse  pressures,  it 
tends  to  predict  high.  (9).   An  interesting  observation 
on  the  use  of  the  Hencky-von  Mises  criterion  is  that  for 
an  unstiffened  cylindrical  shel'l,  Ot  =  2ap.   This  results 
in  failure  of  such  a  shell  occurring  at  aT  =  1.16cfy; 
stated  in  Windenburg's  variables,^  =  1.16. 

The  use  of  formulas  #92  and  #92A  in  the  development 
of  Wenk's  formula  appears  to  be  the  last  time  these 
long-standing  elastic  stress  determinators  were  used  in 
a  yield  criterion.   The  elastic  stress  formulas  of 
Salerno  and  Pulos  (17) ,  having  been  recognized  as  more 
rigorous  and  slightly  more  accurate,  were  used  in  the  next 
collapse  criterion  developed  by  Lunchick.   Lunchick 
utilized  the  Hencky-von  Mises  yield  criterion,  however, 
he  specified  that  collapse  did  not  occur  until  the  load- 
carrying  capacity  of  the  shell  was  completely  exhausted 
at  the  critical  mid-bay  point.   Lunchick' s  three-hinge 
mechanism  stemmed  from  his  assumption  of  plastic  hinges 
forming  at  the  frame  at  the  same  time  a  completely  plastic 
hinge  had  developed  at  mid-bay.   Lunchick ' s  method  predicted 
collapse  pressures  with  better  accuracy  than  any  other 
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in  the  range  of  low  X  values.   However,  its  inclusion  of 
design  curves  to  determine  cylinder  plastic  reserve 
strength  has  made  it  more  difficult  to  apply  in  computer- 
based  design  methods. 

It  is  now  generally  accepted  that  stiffened  cylindrical 
shells  fail    by  one  of  three  distinct  modes  or  any 
combination  of  them.   These  modes  are  (14): 

Mode  1.   Axisymmetric  collapse  of  the  shell  between 

adjacent  ring  frames. 
Mode  2.   Asymmetric  collapse  of  the  shell  between 

adjacent  ring  frames  (lobar) . 
Mode  3.   Overall  asymmetric  collapse  of  the  shell 

and  frames  together. 

Mode  one  had  traditionally  been  referred  to  as 
axisymmetric  shell  yielding  (9,22,23).   It  has  been 
pointed  out,  however,  that  this  mode  is  the  result  of  a 
combination  of  yielding  and  axisymmetric  buckling  (14,22) 
and  the  main  difference  between  this  mode  and  mode  two  is 
the  difference  in  resulting  collapse  patterns,  not  basic 
collapse  mechanisms.   Axisymmetric  collapse  pressures 
have  been  calculated  with  good  accuracy  through  use  of 
yield  criterion  alone  though,  and  it  is  recognized  that 
pre-collapse  stresses,  as  determined  by  von  Sanden  and 


-53- 


Gunter  (20)  or  Salerno  and  Pulos  (17) ,  are  the  primary 
cause  of  failure  in  this  mode. 

The  problem  of  axisymmetric  collapse  by  elastic 
buckling  has  been  solved  by  Salerno  and  Pulos  in  (13) . 
Lunchick  (10)  has  investigated  axisymmetric  failure  by 
inelastic  buckling  and  derived  the  following  expression 
for  it: 
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von  Mises'  and  Windenburg ' s  formulas  for  the  calcula- 
tion of  failure  pressure  by  elastic  asymmetric  (lobar) 
buckling  between  ring  frames  have  been  given.   A  more  recent 
solution  to  this  problem  has  been  found  by  Reynolds  and 
is  outlined  in  reference  (14) .   Reynolds  has  also  developed 
a  formula  for  the  inelastic  buckling  of  stiffened  cylindri- 
cal shells  which  appears  to  be  reasonably  accurate,  but 
is  consistently  high  in  its  calculation  of  collapse  pressures 
(1). 


Mode  thr  e,  collapse  by  general  instability,  has 
been  found  to  be  of  primary  improtance  in  the  preliminary 
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design  of  pressure  vessels  (14).   Solution  of  elastic 
general  instability  collapse  was  first  undertaken  by 
Tokugawa,  who  was  also  the  first  to  identify  this  mode. 
His  formula  for  collapse  pressure,  (which  was  later 
confirmed  by  Bryant),  is  as  follows: 

Eh 


Per 


*m 


(n2  -  1  +  |*)  (n2  +  d2)2 


2    -,  x  Ele 


(n"  -  1) 


R3Lf 


where  X      = 


Lb 


The  original  Tokugawa  formula  as  given  in  (19)  has  a 
number  of  extra  terms  in  it  which  are  seen  to  be 
negligible  for  geometries  of  general  interest. 

The  lengthy  formulas  developed  by  Kendrick  for 
elastic  general  instability  are  given  in  (14),  however, 
a  graphical  solution  of  Kendrick 's  equation  as  developed 
by  Reynolds,  and  refined  by  Ball,  is  included  as  figures 
A-2  and  A-3.   In  this  solution,  the  sum  of  two  effects, 
each  represented  by  a  graph,  is  used  to  determine  collapse 
pressure.   Figure  2  is  the  shell  effect  and  figure  3 
the  frame  effect  which  reflects  the  bending  stiffness  of 
the  cylinder.   Figures  2  and  3  also  indicate  the  Bryant- 
Tokugawa  solution.   Collapse  pressure  by  either  set  of 
curves  is  found  in  the  following  manner: 
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In  order  to  round  out  this  development  of  stiffened 
cylindrical  pressure  vessel  collapse  theories,  mention 
must  be  made  of  the  effect  of  actual  physical  manufacturing 
conditions.   In  all  but  the  smallest  and  most  painstakingly 
fabricated  cylinders,  imperfect  production  processes  are 
the  cause  of  a  number  of  collapse  strength  modifiers. 
Those  modifiers  introduced  into  the  material  itself  are 
in  the  form  of  stresses  due  to  fabrication,  rolling,  and 
welding  of  flat  plating.   In  thick  homogeneous  plating, 
the  homogeneity  itself  may  be  questionable  and  thus  cause 
stress  concentrations.   A  major  strength  modifier  is  the 
imperfect  circularity  of  most  welded  cylinders.   This 
modifier  has  received  some  theoretical  consideration  as 
it  has  a  considerable  effect  on  the  resistance  of  a  stiff- 
ened cylinder  to  failure  by  general  instability.   A 


formula  for  this  is  due  in  part  to  Kendrick  (2) 

85 


P1R         1  +  ~B —     Ec     o           /  Pl 
b  l-LJL ]  +  ££_  (n2-l)e/ 

Ab  1+0  Rl\  lPcrPr 


(—) 

I  P„T.p,/ 


As  long  as  the  a max  is  below  material  yield  stress, 
collapse  by  general  instability  is  not  indicated. 
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This  section  has  been  included  to  give  a  brief 
background  of  the  development  of  theoretical  collapse 
predictors  for  stiffened  cylindrical  shells.   In  particu- 
lar, a  history  of  the  variables  of  X   and  ty     is  helpful 
in  understanding  their  use  in  this  paper. 
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APPENDIX  B 
DISCUSSION  OF  VERY  LOW  LAMBDA  DESIGN 

As  noted  in  Section  IV,  overall  structural  efficiency 
approaches  its  highest  value  as  \p  goes  to  1.0  and  A  goes 
to  0.   For  T-shaped  frames  of  the  geometry  used  in  this 
report,  X  values  in  the  area  below  .2  may  be  considered 
very  low  lambda  design.   For  other  frame  geometries,  this 
may  be  higher.   For  example,  the  sandwich  construction 
in  figure  A-4  has  a  X  value  of  about  .337.   The  frames  used 
there  have  thicker  webs  than  flanges,  whereas  the  reverse 
is  true  for  the  frames  of  this  report.   The  membrane  shell 
design,  and  other  complications  resulting  from  the  use  of 
high  strength  titanium  also  plays  a  part  in  deciding  the 
optimum  X  value. 

It  was  noted  in  Section  IV  that  the  danger  of  general 
instability  collapse  would  be  greatly  increased  by  the 
use  of  a  very  low  lambda  design.   This  has  been  shown  to 
be  true  in  model  tests  of  cylindrical  shells  with 
geometries  similar  to  that  of  figure  A-4.   In  fact,  in  the 
testing  of  these  models,  a  new  mode  of  collapse  was  observed 
-  that  of  plastic  general  instability  -  and  new  theories 
had  to  be  developed  to  predict  it.   It  may  also  be  observed 
that  collapse  of  these  models  took  place  at  ty   value  in 
excess  of  3.0.   This  suggests  that  at  very  low  lambda 
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values,  <C  i    becomes  a  rather  meaningless  term  unless 
the  framing  is  very  small  in  comparison  to  the  shell 
thickness.   Indeed,  the  large  effect  of  framing  at  these 
A  values  throws  into  doubt  the  meaningf ulness  of  A  itself, 
This  doubt  has  been  substantiated  but  not  fully  explained 
in  reference  (6) . 
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WELDABLE   TITANIUM  JACKET   <J     =    120,000   psi 


TITANIUM  ALLOY 
gy  =   175,000  psi 
(unweldable) 


SCHEMATIC   SKETCH    of   an  OCEANOGRAFHIC    VEHICLE 
WITH     •    9-FOOT  OUTER   DIAMETER   from  DTMB    1677 


FIGURE     A -4 
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APPENDIX  C 

FRAME  PARAMETERS 

A  T-shaped  frame  was  used  in  the  calculations 
leading  to  figures  5-10.   The  proportions  of  this 
frame  were  taken  from  a  work  by  McGinley  (11)  and  are 
as  follows: 

frame  width/frame  depth  =  1.0 
frame  width/web  thickness  =  8.5 
flange  thickness/web  thickness   =   1.7 

Figure  A— 5  is  a  graph  of  b/D  versus  X  for  the  frames 
for  both  formulas  #92A  and  Wenk's.   Since  the  geometry 
of  optimized  vessels  was  found  to  be  a  function  of  required 
collapse  pressure, and  not  required  diameter,  the  optimized 
frame  dimensions  should  be  proportional  to  the  major 
cylinder  dimensions.,  hence  the  b/D  parameter.   Figure  A-4 
may  be  used  in  conjunction  with  figures  8,  9,  and  10  to 
determine  the  frame  dimensions.   Care  should  be  taken, 
however,  to  match  yield  formulas  and  frame  factors  when 
using  these  figures  together. 


.360 


J.35G-. 


\ 


lOOb/D    vs 
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-64- 

APPENDIX  D 
EXAMPLES  OF  DESIGN  ITERATION 
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Example  2 : 

This  example  is  included  to  indicate  the  results 
of  the  iterative  method  when  applied  to  relatively  low 
pressure  vessels. 

Given:       Use  of  #92A   Frame  Factor  5 
cy   =   30,000  psi 
pc   =   400  psi 
D    =   40  ft. 


Steps: 


1.  Initial  h/Dm  guess  from  hoop  stress 
(h/Dm)   =   400/(2*30,000) 

(h/Dm)   =   .00666 

A        =   .495 

ij>       =   1.30 

pc      =   1.30*400   =   520  psi 

which  is  too  high 

2.  Reduce  h/Dm  by  division 

pcl     =    400      =   308  psi 

1.3 

(h/Dm)   =   .00512 

X        =   .565 

ip  =   1.39 

p        =   1.39*308   =   428  psi 

which  is  still  too 

high 

*act    =    428      =   1'07 
400 


Answer: 


Results: 
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Reduce  h/Dm  again 

p^9      =    308      =   288  psi 
°  "TToT 

(h/Dm)   =   .00480 

X       =   .584 

ip  =   1.41 

pc      =   1.41*288   =   406  psi 

which  is  close  enough 


■X 

= 

.584 

h/D 

= 

.00480 

f\ 

= 

1.27*30,000 
80,000 

n 

= 

.476 

* 

= 

1.41 

D 

= 

40  ft. 

h 

= 

2.3  in. 

L 

= 

4.5  ft. 

Wt 

= 

Pc 

A 

3720n 

Wt 

= 

.23 

406 


3720(.476) 
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Example  3: 

This  example  is  exactly  like  example  2  except  that 
a  higher  strength  material  is  used. 

Given:        Use  of  #92A   Frame  Factor  5 

ay  =  50,000  psi 

pc  =  400  psi 

D  =  40  ft. 


Answer: 

X 

= 

.825 

h/°m 

= 

.00232 

n 

= 

1.11*  50,000 
80,000 

n 

= 

.694 

* 

= 

1.735 

Results: 

D 

= 

40  ft. 

h 

= 

1.11  in. 

L. 

= 

1.83  ft. 

Wt 

= 

.157 

A 
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Example  4: 

This  example  is  included  to  indicate  the  use  of  the 
frame  parameter  curves  in  Appendix  C. 

Given:        Use  of  #92A   Frame  Factor  6 

a  =  100,000  psi 

D  =  40  ft. 

pc  =  1500  psi 


Steps: 


1.  Initial  h/Dm  guess 

(h/Dm)    =   1500/(2*100,000) 
(h/Dm)    =   .007 

X  =   .48 

if>       =   1.33 

pc       =   1.33*1500   =   2000  psi 

2 .  Reduce  h/Dm 
1500 


Pel 

1.33       =   ±1JU 

(h/Dm) 

= 

1130/(2*100,000) 

(h/Dm) 

= 

.0056 

X 

= 

.54 

i> 

= 

1.40 

Pc 

= 

1.40*1130   =   1580  psi 

*act 

= 

1580        _ 

1500           -L,U:3 
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3.   Reduce  h/Dm- again 


pc2   =  T705       =   1075  psl 
(h/Dm)  =   1075/(2*100,000) 


(h/Dm)  = 

.00537 

X 

.55 

1i 

1.41 

Pc 

1.41*1075 

=    1515  psi 

- 

close  enough' 

Answer:        X      =   .55 


(h/Dm) 

— 

.00537 

n 

= 

100,000 
80,000 

* 

= 

1.41 

b 
D 

= 

.00248 

Results:       D 

= 

40  ft. 

h 

= 

2.6  £n. 

L 

= 

1.42  ft. 

Wt 

A 

= 

.252 

Frame  Dimensions: 

b 

= 

1.18  in. 

w 

= 

8.05  in. 

f 

= 

10.1  in. 

t$ 

= 

2.05  in. 

1.59 
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APPENDIX  E 
FORMULAS  AND  DERIVATIONS 

In  this  appendix,  the  formulas  used  in  the  computer 
programs  in  Appendix  F  are  listed.   Also,  the  derivation 
of  formulas  (A)  and  (B) ,  Section  II,  and  the  solution 
of  the  Hencky-von  Mises  criterion  for  unstiffened  tubes 
are  presented. 

FRAME  FORMULAS:   The  formulas  below  are  used  to 
predict  frame  failure  by  instability.   Collapse  of  the 
frame  by  yielding  was  not  investigated  as  the  frames 
were  not  required  to  hold  up  after  shell  failure  and  the 
maximum  stress  in  the  frame  is  always  less  than  the 
maximum  stress  in  the  shell  (19)  .   However,  this  stress 
relationship  may  not  hold  in  the  very  low  lambda  region. 


Foppl  Formula: 

24EIb 

Pc   -   ; 

Db3L! 

Formula  #8  8 

von  Sanden  and  Gunther: 

24EIb  (1+6) 
P, 


c 


Db3b  (1  +  ^|i-) 


u     n      bh 
where  B  ■  

A  +  bh 
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L/Dm 

e  =  i.82   h/Dm 

3=2^  NB 


■b   9 

N  =  cosh9  -  sin( 


sinhB  +  sin0 


SHELL  FORMULAS:  The  first  formula  below  predicts 
collapse  of  the  shell  between  frames  due  to  instability, 
The  next  four  predict  collapse  due  to  yielding  of  the 
material  according  to  different  yield  criteria. 

Formula  #92 

von  Sanden  and  Gunther 

Failure  due  to  axial  stress  exceeding  yield  stress: 

^ha/        1       \ 


K  = 


Dm  Y   I    .5  +  1.81KC 

sinh8  -  sin8 
sinh0  +  sin8 

r   -  -85  "  B 

i  +  e 

B  and  B  are  defined  in  Formula  #88 

Formula  #92A 

von  Sanden  and  Gunther 

Failure  due  to  circumferential  stress  exceeding  yield 

stress: 

h     /      1      \ 


pc    ■    V 


1  +  HC 
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(1  +  M)  sinhy  cosy  +  (,1-M)  coshy  siriy 
H  =  -2 - — 

sinh6  +  sin9 


M-   ^L 


1  -  P2 


C  is  defined  in  Formula  #92 

Wenk's  Formula 

This  formula  was  actually  developed  by  Wenk  and  others 
at  DTMB  using  #92  and  #92A  to  determine  the  principle 
stresses.   These  were  applied  at  mid-bay  and  mid-thickness 
using  the  Hencky-von  Mises  yield  criterion. 


Pc       DT°y 
m 


F2  -  ,5F  +  .25 


F  =  1  +  2QC 

.  ,  e    e     ,e.e 

sinhw-  cosj  +  coshj   sin~- 
sinh0  +  sine 

6  and  C  are  defined  in  Formula  #92 

Lunchick's  Formula 

This  formula  uses  the  principle  stress  relations  developed 
by  Salerno  and  Pulos  rather  than  #92  and  #92A  (9).   It 
applies  them  to  the  mid-bay  location  using  the  Hencky- 
von  Mises  criterion,  however,  failure  is  not  predicted 
until  the  yielding  which  occurs  first  in  the  outer  surface 
develops  into  a  fully  plastic  hinge.   Lunchick's  formula 
predicts  this  yield  pressure,  but  curves  (included  here 
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as  figure  A-6)  must  be  used- to  determine  the  actual 
collapse  pressure. 

1 


o  h 


X2  -  xz  +  z2 


X  =  1  +  HC 

Z  =  .5  +  3.62TC 

,  e    e     .e  .  e 

sxnhT  cos?  -  cosh7  sin7 

m  Si  _  £ £ £ ^ 

sinh9  +  sinB 
6,  C,  and  H  are  defined  in  Formula  #92"A 


Entry  into  the  curves  requires  computation  of  the 

following  ratios: 

Circumferential  Bending  Stress      qBT 
Circumferential  Membrane  Stress     o^, 

Lonaitudinal  Membrane  Stress       aF 


Circumferential  Membrane  Stress     T 


aF 


Dm  Pc 
4h 


aT  =   2aF(l  +  2QC) 


aBT  =   2aF  TC 


•74- 


d  5s 


<D 

^i 

D 

\n 

<D 

w 

^ 

O 

D 

u 

n 

~m 

w 

<u 

1) 

u 

i/) 

•A, 

^c 

cc 

■a 

^ 

^-< 

r-( 

<D 

0 

■H 

u 

>• 

1.28 


1.24 


1.20 


1.16 


1.12 


1.08 


1.04 


1.00 


Circumferential   Bending   Stress 
Circumferential  Membrane  Stress 


Curves    for  Lunchick's   Collapse   Pressure 
Factor   -   From  DTMB  Report    1291 


FIGURE     A-6 
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GENERAL  INSTABILITY  FORMULAS:   While  Tokugawa ' s 
formula  and  graphs  of  Kendrick's  third  solution  for 
general  instability  have  been  given  in  Appendix  A, 
Eryant ' s  simplification  of  Kendrick's  first  solution 
has  been  used  in  the  computer  program. 


Bryant's  Formula: 


cr 


2JLe 


X* 


(n2  -  1  +   |  ) (n2  +  X2) 


n2  -  1 


-EI- 


R  aL, 
m   1 


ttR 


an  integer 


n  =   the  integer  which  gives  the  lowest  pcr. 


Since  out  of  roundness  effects  general  instability 
collapse,  a  formula  taken  from  (2)  has  been  used  in  the 
following  form: 


-max 


P  R 

^c 


1  + 


85g 
B 


1+6 


Ec 
R,  : 


'(n2  -l)e 


'cr 


"  P, 
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e  =  eccentricity 

p  =  desired  collapse  pressure 

n,  pcr,   3/  and  B  are  determined  in  previous 
formulas . 
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X       FROM  PRINCIPLE  STRESS  CRITERION  AT  MID-BAY  EXTERIOR, 

m  ' — ■ — ■ 


FORMULA  #92A 


1.   Formula  #92A 


2.      ty 


2  hgy 

Dm 


1    +    H    I    -85    -    B 


V       1  +  3 


o    h   g 


T|>  1     @        1  *        1 

'.8  5    -    B    \ 


1    +    H 


1   +   6 


1    +    H   (  ,85    "    B     1 
V"  1  +   3 


H     f-85    "    B     )  =       0 

V      1    +    3    / 

H  =    0 

3    sinh^-  cos^  +   coshS-  sinS- 
H   =    —      


sinh      +    sin 
H=0@6      =3.74    radians 
9      =       (1.82    L/Dm)  /  \|h/Dm 


L/Dm=      2.055  \JT7d; 
h/Dm=  .23  6(L/Dm): 
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3.       A      = 


^m  - 


\ 


(L/Dn) 2          i 

°Y 

(h/Dm) 3           | 

E 

2.055    Dm 

°y 

h 

E 

Equation    (A) ,    Section   II 
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Xm  FROM  HENCKY-VON  MISE5  CRITERION  AT  MID-BAY  AND  MID- 


THICKNESS  , 

VON  SANDEN  AND  GUNTHER  FORMULAS 


Wenk' s  Formula 
2hav 


D 


m 


IF*    -    F    +    1 
»  2         4 


1    +    2Q 


/  .85    -    B   \ 
\        1    +    6  / 


£         i  £         i 

sinh   2    cos    2    +   cosh   2    sin   2 
sinh   6    +    sin   8 


2.       \\i 


* 


2A  av 

1 

■.= 

1-154@   \L2 

2         4 

1.154 

F2 

-  n+  i 

2         4 

= 

.75 

F      =    ,25±    .75 


.*.    1      =    1    +    2Q 


35    -    B 


1    +    6   / 
.     Q      =    0    @    6    =    4.68    radians 


(1.82    L/Dm)    /  \[h/Dm 
L/Dm   =    2.565    Jh/Dm 
h/D^    =       ,152     (L/Dm): 
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3.  X 


(L/Dm) 


(h/Dm) 


E 


m 


\ 


2.565  D 


m 


E 


Equation  (B) ,  Section  II 
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SHELI  IDELY  SPACED  STIFFENERS 

1.   He    --von  Mises  Criterion 

C;l       :'2      ~   three  principle  stresses 

2c        (a i   -a2    )2  +  (a2  -a3  )2  +  (a3  -"C^)2 


For 

1 

anstiffened  tube: 

al 

= 

2a2 

a3 

2c, 

= 

a22  +  a22  +  4a22 

oy- 

= 

3a| 

er 

- 

1 

9 

1.732  a 

a 
1 

2 

1.732  ay 

.577ay 


1.154ay 
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APPENDIX  F 
COMPUTER  PROGRAMS 

Two  computer  programs  were  used  to  produce  theoretical 
data  for  this  report.   Program  1  utilizes  equation  (A) 
of  Section  II  to  determine  curves  of  h/Dm  versus  Xm 
for  eight  values  of  Oy.   Since  equation  (A)  is  based  on 
the  9  value  at  which  the  variable  H  in  formula  #92A 
goes  to  zero,  the  program  first  solves  for  that  6  and 
then  gees  on  to  vary  h/Dm  and  Cy,    solving  for  Xm.   This 
program  was  used  only  to  define  the  upper  limits  of  the 
yield  failure  area  and  so,  is  not  actually  necessary  for 
the  iterative  optimization  proceedure  of  Section  IV.   A 
flow  chart  and  variable  chart  are  included  with  the  program 
listing  and  printout. 

Program  2  follows  a  scatter-gun  philosophy  in  that 
it  determines  collapse  pressures  of  given  geometries  for 
all  three  modes  of  failure  by  various  formulas.   All  of 
the  formulas  are  mentioned  in  Appendix  A  and  completely 
defined  in  Appendix  E.   The  input  of  this  program  consists 
of  outside  diameter,  desired  range  of  collapse  pressures, 
frame  proportions,  material  properties,  and  range  of 
values.   Shell  thickness  is  determined  by  the  hoop  stress 
relation  and  frame  size  by  a  modified  Foppl  formula.   Program 
2  was  not  intended  to  be  used  to  produce  a  complete 
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optimized  design  for  a  given  set  of  requirements.   Its 
main  use,  in  solving  for  and  printing  collapse  pressures 
for  all  modes,  is  pointing  up  trends  and  areas  where 
various  formulas  fail.   It  has  provisions  for  any 
number  of  solutions  to  be  printed  and  any  geometry  to 
be  investigated.   In  short,  it  is  meant  to  be  used  as 
a  tool  by  a  designer  who  optimizes,  using  judgment  and 
experience  as  well  as  computer  time.   A  flow  chart  and 
variable  chart  are  included  with  the  program  listing 
and  a  sample  of  its  printout. 
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FROG  RAM   1 

THIS    FROGRAM      DETERMINES    LAMBDA    MAX    FOR  CHANGING   H/D 
AND    YIELD   STRESS    USING  A    FORMULA    DERIVED    FROM   92A . 

P2   =    .09 
NUM  =    1 
NEND    =    50 
TH  =   3.7 
DTHP   =    .1 
DTHM  =    .1 
1    NUM  =   NUM  +    1 

IF(NUM.EO.NEND)   GO  TO   30 
THT  =   TH/2. 

DEN    =   SINH(TH)+SIN(TH) 

EHO  =    (l.+(3*P2/(l-P2))**.5)*SINH(THT)*COS(THT). 
EH1    =    (l.-(3*P2/(l-P2))**.5)*COSH(THT)*SIN(THT) 
EH   =    -2*(EH0+EH1)/DEN 
IF(EH)    19,30,9 
9  NP   =   NF+1 

IF(NF.EQ.3)   GO   TO    10 
DTHP   =   DTHP/2 
10  TH   =   TH-DTHM 
M  =    1 
GO   TO    1 

19  M  =  M+l 
IF(M.EQ.3)    GO   TO   20 
DTHM  =   DTHM/2 

20  TH   =   TH   +   DTHP 
NP  =    1 

GO      TO    1 
30  WRITE(6,180)   NUM,EH,TH 
180  FORMAT (»    AFTER',  13,  •    ITERATIONS,    H    IS', F8.6    •    AND   THETA 

E  =   30000000.0  IS'    F6  3) 

DOH  =   0.0  *       * 

DDOH   =40.0 

NDOH   =15 

DY  =   10000.0 

NY  =   8 

DO    50  N    =    1-NDOH 

WRITE(6,205)   DOH 

Y  =  30000.0 

DO  40  M  =   1,NY 

ELAM  =    (TH/1.82*DOH*Y/E)**.5 

WRITE(6,210)   Y,ELAM 
40  Y   =   Y+DY 
50  DOH   =   DOH.+    DDOH 
205   FORMAT (•    LA  BDAS    FOR  DM/H   =    'F7.2    '    ARE    •) 
210  FORKAT(2Fl2.4) 

END 
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FRINTOUT  OF    PROGRAM   1       (Sample) 


AFTER    50    ITERATIONS 
LAMBDAS    FOR   DM/H   = 


H    IS -.000000  AND   THETA    IS    3.74 
0.0     ARE 


30000.0000 

0.0 

40000.0000 

0.0 

50000.0000 

0.0 

60000.0000 

0.0 

70000.0000 

0.0 

80000.0000 

0.0 

90000.0000 

0.0 

100000.0000 

0.0 

LAMBDAS  FOR 

DM/H 

40.00  ARE 

30000.0000 

0.2867 

40000.0000 

0.3311 

50000.0000 

0.3702 

60000.0000 

0.4055 

70000.0000 

0.4380 

80000.0000 

0.4683 

90000.0000 

0.4967 

100000.0000 

0.5235 

LAMBDAS  FOR 

DM/H 

80.00  ARE 

30000.0000 

0.4055 

40000.0000 

0.4683 

50000.0000 

0.5235 

60000. 0000 

0.5735 

70000.0000 

0.6194 

80000.0000 

0.6622 

90000.0000 

0.7024 

100000.0000 

0.7404 

LAMBDAS  FOR 

DM/H 

=   120.00  ARE 

30000.0000 

0.4967 

40000.0000 

0.5735 

50000.0000 

0.6412 

60000.0000 

0.7024 

70000.0000 

0.7587 

80000.0000 

0.8110 

90000.0000 

0.8602 

100000.0000 

0.9068 

LAMBDAS  FOR 

DM/H 

=   160.00  ARE 

30000.0000 

0.5735 

40000.0000 

0.6622 

50000.0000 

0.7404 

60000.0000 

0.8110 

70000.0000 

0.8760 

80000.0000 

0.9355 

90000.0000 

0.9933 

100000.0000 

1.0470 
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(  START  J 


P2=.09  TH=3.7 
NUM=1  NEND=50 
DTHP=DTHM=.l 


CALCULATE 
EH  (H) 


NP  =  NP  +  1  I 


is 


YES 


WRITE  NUM 
EH  TH 


E,  DOH,  DDOH 
NDOH,  DY,  NY 


\^^  M>2  ? ~^ 

NO 

DTHM 

dp 

DTHM 

2 

TH=TH+DTHP 

TH=TH-DTHM  i 
/ 


J 


/ 


Y  =  30,000 


/ 


i  CALCULATE  ELAM. 
^WRITE  Y,  ELAM 


FLOW  CHART 
PROGRAM  1 
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NOMENCIATURE  OF  PROGRAM  1 


DOH 

DDOH 

NDOH 

E 

EH 

ELAM 

P2 

TH 

Y 

DY 

NY 


m 
Incremental  change  in  Dm/h 

Number  of  changes  in  Dm/h 

Young's  Modulus 

Variable  H  of  yield  failure  formulas 

*m 

Poisson's  Ratio  squared 

Variable  8  of  yield  failure  formulas 

Material  ay 

Incremental  cnange  in  Oy 

Number  of  changes  in  Gy 
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FROGRAM  2 


THIS    PROGRAM  DETERMINES    THEORETICAL  COLLAPSE    PRESSURES 
FOR  STIFFENED   CYLINDRICAL    PRESSURE   VESSELS.    THREE   MODES 
OF  COLLAPSE  ARE    INVESTIGATED  WITH    EMPHASIS    ON    YIELD 
PRECIPITATED    BUCKLING    BETWEEN    THE   FRAMES.    FRAME   SIZE 
IS    DETERMINED    USING  A    MODIFIED    FOPFL  FORMULA.    ALL    INPUTS 
ARE  HANDLED   BY  MAIN    PROGRAM  CARDS    -  NO  DATA   CARDS   ARE 
USED. 

MF  =   6 

D  =  480.0 

PC   =   300.0 

DO   700  NU  =   1,14 

Y  =   80000.0 

E  =   30000000.0 

Rl   =   1.0 

R2   =    .11764 

R3  =   1.7 

P  =    .3 

F2   =  P**2 

SR  =    .38 

IF(FC.LT. 1000.0)   SR  =    .46 

DSR   =    .02 

NSR   =   20 

T   =   D*FC/(2*Y) 

DM  =   D-T 

TOD   =   T/DM 

WRITE (6, 160)    T,TOD 

DO   500   IN   =   1,NSR 

EL  =   SR**2*DM*E/Y*(T/DM)**I.5 

FRAMD1    =    12.0 

DFRP   =1.0 

DFRM  =1.0 

NUM  =    1 

NEND   =    70 

RATIO   =   EL*FC/(24*E*MF) 

NUM  =   NUM  +    1 

IF(NUM.EQ.NEND)   GO  TO   30 

F  =   FRAMDl*Rl 

B  =   FRAMDl*R2 

T2   =   B*R3 

W   =   FRAMDl    -  T2 

YB   =   YENTRD(B,T,B,W,T2,F) 

YBI   =   YMOM(B,T,B.W,T2,F>YB) 

DBT  =    (D-2*(YB+T))**3 

RAT   =   YBI /DBT 

IF(RAT-RATI,)    19,30,9 

NP   =   NP+1 
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IF(NP.EQ.3)   GO   TO    10 
DFRF   =   DFRF/2 
10  FRAMDl    =   FRAMDl    -  DFRM 
M  =   1 
GO   TO    1 

19  M  =   M+l 
IF(M.EQ.3)   GO  TO  20 
DFRM  =  DFRM/2 

20  FRAMDl    =   FRAMDl   +   DFRP 
NP  s    1 

GO   TO   1 
30  YF   =   YENTRD(0.0,0.0,B,W,T2,F) 
Yl    =   YENTRD(EL,T,B,W,T2,F) 
YI   =   YMOM(0.0,0.0,B,W,T2,F,YF) 
Y1I   =  YM0M(EL,T,B,W,T2,F,Y1) 
YB10YI   =   YBI/YI 
Y1I0YI   =   Y1I/Y1 
RONE   =    (D-2*(T+Yl))/2 
DB   =   D-2*(T+YB) 
A   =   F*T2   +   B*W 
AB   =  A+B*T 
EB  =   B*T/AB 
THETA  =   10*(12-12*P2)**.25 

TH   =   THETA*EL/DM/(100*T/DM)**.5 

THT  =  TH/2 

DEN   =   SINH(TH)+SIN(TH) 

EN   =    (COSH(TH)-COS(TH))/DEN 

EK   =    (SINH(TH)-SIN(TH))/DEN 

BETA   =   2.*EL*EN*EB/(B*TH) 

EHO  =    (l+(3*F2/(l-F2))**.5)*SINH(THT)*COS(THT) 

EH1    =    (1-(3*F2/(1-F2))**.5)*C0SH(THT)*SIN(THT) 
EH   =    -2*(EH0+EHl)/DEN 

EQ   =    -(SINH(THT)*COS(THT)+COSH(THT)*SIN(THT))/DEN 
ET   =   -(SINH(THT)''C0S(THT)-G0SH(THT)*S1N(THT))/DEN 
SUBDEN   =    (.85-EB)/(l+BETA) 
EF   =   1.    +    2*EQ*SUBDEN 
EK  a   1.    +EH*SUBDEN 
EZ   =    ,5+3.62*ET*SUBDEN 
PTOK  =   2  4*E*Y1 1 / ( DM**3*EL* ( 1 - P2 ) ) 
PWIND   =  2 ,6*E*(T/DM)**2 ,5/(EL/DM-.45*(T/DM)**.5) 
DENONE   =    (1+(.85*BETA/EB))/(1+BETA) 
DENUM  =   2*T*Y/DM 

FY88   a   24*E*YBl/(DB**3*B*DENONE) 
•     SI92    =    l/(.5+1.81*EK*SUBDEN) 
S192A    =    l./(l+EH*SUBDEN) 
SIWENK  =   l./(EF**2-EF/2+.25)**.5 
F  LUNCH  =   2*T*Y/D / ( EK**2 -EX*EZ+EZ**2 )** . 5 
FBRY  =  2*E*T/DM/850.   +   64*E*Y1I/(DM**3*EL) 
SIGOOR   =    PC  --D*B*DEi;ONE/(AB*2) 

S1GOOR   =   SIGOOR  +    E*( (D-DB)/2 . )*4*(IC/(F BRY-FC) )/RONE**2 
S1GF   =    1000*DM/(4*T) 
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SIGT  =  SIGF*2*(1.+2*EQ*SUBDEN) 

SIGBT   a   SIGF*2*ET*SUBDEN 

SIGFOT   =   SIGF/SIGT 

SIGBOT   =   SIGBT/SIGT 

R   =   D/2. 

FI   =   3.14159 

VT  =   FI*(EL*(R**2-(R-T)**2+B*((R-T)**2-(R-T-W)**2) 
1      +F*(  (R-T-W)**2-(R-T-W-T2  )**2)) 

WT   =   VT*. 284*12 /EL 

VOL  =   PI*(R/12.)**2*12 

SF92A    =   SI92A*DENUM*V0L/(WT*1 00000) 

S  FWENK   =  S IWENK*DEKUM*VOL/ (WT*1 00000 ) 

Write   statements    for  whatever  solutions   are  required 
500  SR   =  SR+CSR 
700   FC   =   FC+100.0 

Format   statements 

END 


FUNCTION  SUBROUTINES 


FUNCTION  YENTRD(EL,T.B,W,T2,F) 

902  TOP  =  2*F*T2*W  +  F*T2**2  +  B*W**2  -  EL*T**2 

903  BOT  =  2*(F*T2  +B*W  +  EL*T) 

904  YENTRD  =  TOP/BOT 
RETURN 

END 


FUNCTION  YMOH(EL ,TLBLWLT2 ,F , Y) 

922  FIRST  =  (EL*T**3  +  B*W**3  +  F*T2**3)/12 

923  SECOND  =  EL*T*(T/2+Y)**2  +  B*W*(Y-W/2 )**2 

924  THIRD  =  F*T2*(T2 /2+W-Y)**2 

925  YMOM  ■  FIRST  +  SECOND  +  THIRD 
RETURN 

END 
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FRINTOUT  OF   PROGRAM  2       (Sample,    continued    on   next    page) 

T    IS      3.0 INCHES   AND   T/DM   IS    0.00529 

LAMBDA, EL, FRAME   FACTOR  0.380         12.883  4 

SECOND   MOMENTS :F,B,1        200.436        439.394         1481.279 

F  AREA         WINDENBG        FORM88         FORM92  FORM92A  WENK 

28.2750   -28193.383        590.321      1.857  1.488  1.665 

LAMBDA, EL, FRAME   FACTOR  0.400        14.275  4 

SECOND  MOMENTS :F,B,1        225.964        485.456        1680.225 

F  AREA         WINDENBG         FORM88         FORM92  FORM92A  WENK 

30.022      -42473.102        609.692      1.830  1.469  1.651 

LAMBDA, EL, FRAME   FACTOR  0.420         15.738  4 

SECOND   MOMENTS:F,B,l        253.112        533.708        1893.137 

F  ARAE        WINDENBG         FORM88         FORM92  FORM92A  WENK 

31.774     -90844.625        628.875      1.801  1.451  1.637 

LAMBDA, EL, FRAME   FACTOR  0.440         17.273  4 

SECOND   MOMENTS :F,B,1         281.881         584.134        2120.124 

F  AREA        WINDENBG         FORM88         FORM92  FORM92A  WENK 

33.531        467240.063     647.914      1.770  1.434  1.624 

LAMBDA, EL, FRAME   FACTOR  0.460        18.879  4 

SECOND  MOMENTS :F,B,1        312.279        636.729        2351.309 

F  AREA        WINDENBG         FORM88         FORM92  FORM92A  WENK 

35.293        62893.883        666.852      1.736  1.420  1.612 

LAMBDA, EL, FRAME  FACTOR  0.480        20.556  4 

SECOND  MOMENTS :F,B,1        344.308        691.478        2616.743 

F  AREA        WINDENBG         FORM88        FORM92  FORM92A  WENK 

37.058        33034.949        685.733      1.700  1.400  1.599 

LAMBDA, EL, FRAME   FACTOR  0.500        22.305  4 

SECOND   MOMENTS :F,B,1         377.968        748.369        2886.512 

F  AREA        WINDENBG        FORM88         FORM92  FORM92A  WENK 

38.828        22097.688        704.602      1.662  1.383  1.587 

LAMBDA, EL, FRAME   FACTOR  0.520        24.125  4 

SECOND   MOM£NTS:F,B,l        413.264        807.395        3170.714 

F  AREA        WINDENBG         FORM88        FORM92  FORM92A  WENK 

40.600        16434.441        723.510     1.622  1.367  1.575 

LAMBDA, EL, FRAME   FACTOR  0.540        26.016  4 

SECOND  MOMENTS :F,B,1        450.198        868.543        3468.377 

F  AREA         WINDENBG        FORM88         FORM92  FORM92A  WENK 

42.376        12977.973        742.510     1.581  1.350  1.563 

LAMBDA, EL, FRAME   FACTOR  0.560        2  7.979  4 

SECOND   MOMENTS ,F,B,1         488.769        931.803        3782.55 

F  AREA         WINDENBG         FORM88        FORM92  FORM92A  WENK 

44.154        10653.945        761.666      1.538  1.333  1.551 
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FLANGE, WEB: LENGTH, THICKNESS         9.80      1.96      7.84  1.15 
LUNCHICK  RATIOS :F/T,B/T                    .750      .0080 

LUNCH ICK  BRYANT  OUT   ROUND        WT/FT        EF92  EFWENK 

1.629  2478.0  54958.6      26118.07    1.246  1.327 

FLANGE, WEB: LENGTH, THICKNESS      10.10      2.02      8.08  1.19 
LUNCHICK  RATIOS, F/T,B/T                     .742       .0095 

LUNCHICK  BRYANT  OUT  ROUND        WT/FT        EF92  EFWENK 

1.610  2526.2  55407.4     25657.38   1.253  1.339 

FLANGE, WEB: LENGTH, THICKNESS      10.39      2.08      8.31  1.22 
LUNCHICK  RATIOS :F/T,B/T                     .735       .0112 

LUNCHICK  BRYANT  OUT  ROUND        WT/FT        EF92  EFWENK 

1.592  2571.9  55852.6     25234.33    1.257  1.348 

FLANGE, WEB: LENGTH, THICKNESS      10.68      2.13      8.54  1.26 
LUNCHICK  RATIOS :F/T,B/T                     .727       .0131 

LUNCHICK  BRYANT  OUT  ROUND        WT/FT         EF92  EFWENK 

1.574  2615.4  56292.1      24844.41    1.265  1.356 

FLANGE, WEB: LENGTH, THICKNESS      10.95      2.19      8.76  1.29 
LUNCHICK  RATIOS :F/T,B/T                     .720      .0151 

LUNCHICK  BRYANT  OUT  ROUND        WT/FT        EF92  EFWENK 

1.556  2656.7  56722.0     24483.86   1.267  1.363 

FLANGE, WEB:  LENGTH, THICKNESS      11.22      2.25      8.98  1.32 
LUNCHICK  RATIOS :F/T,B/T                     .713       .0174 

LUNCHICK  BRYANT  OUT  ROUND        WT/FT        EF92  EFWENK 

1.538  2695.9  57140.1      24149.51    1.268  1.369 

FLANGE, WEB: LENGTH, THICKNESS      11.49      2.30      9.19  1.35 
LUNCHICK  RATIOS :F/T,B/T                     .706      .0198 

LUNCHICK  BRYANT  OUT  ROUND        WT/FT        EF92  EFWENK 

1.519  2733.4  57544.6     23838.47    1.269  1.374 

FLANGE, WEB: LENGTH, THICKNESS      11.75      2.35      9.40  1.38 
LUNCHICK  RATIOS :F/T,B/T                    .700      .0225 

LUNCHICK  BRYANT  OUT  ROUND        WT/FT        EF92  EFWENK 

1.501  2769.0  57932.6     23548.48    1.269  1.378 

FLANGE, WEB:  LENGTH,  THICKNESS      12.00      2.40      9.60  1.41 
LUNCHICK  RATIOS :F/T,B/T                     .694       .0253 

LUNCHICK  BRYANT  OUT  ROUND        WT/FT        EF92  EFWENK 

1.482  2803.1  58301.6     23277.38   1.267  1.379 

FLANGE, WEB: LENGTH, THICKNESS      12.25      2.45      9.80  1.44 
LUNCHICK  RATIOS :F/T,B/T                     .687       .0284 

LUNCHICK  BRYANT  OUT   ROUND        WT/FT         EF92  EFWENK 

1.462  2835.6  58649.1      23023.39   1.265  1.381 
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(START         ") 

MF,    D,     PC 

—    - 

Y,    E,     Rl,     R2,     R3,    P, 
SR 

CALCULATE  T,  DM,  TOD 
WRITE  T,  TOD 


r 


-\ 


CALCULATE  EL,  RATIO 


FRAME1,  DFRP,  DFRM,  NUM, 
NEND 


1A 


is 
NUM=NEND 


NO 


YES 


CALCULATE  F,  B,  T2 ,  W, 
YB,  YBI,  DBT,  RAT 


is 

YT^RATIO 

NO 

is  " 
RAT>RATI0 


NO 


YES 


2A 


FLOW  CHART 
PROGRAM  2 
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All  calculations  used  to 
find  solutions  to  various 
formulas  for  pc  in  differ- 
ent modes 

- 

All  necessary 
write  statements 


© 


END 


J 


1  — I 


FRAMD1=F RAMD 1 - DF  RM 


L-jyJ 


1A 


Loop  2  to  2A  is  the  same  as  1  to  1A  in  its  method 
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NOMENCLATURE  OF  PROGRAM  2 


B      -  Width  of  frame  in  contact  with  shell 

D      -  Outside  diameter  of  shell 

DSR    -  Incremental  change  in  slenderness  ratio 

E      -  Young's  Modulus 

EL     -  Length  of  shell  between  frames 

F      -  Frame  overall  width 

FRAMDl  -  Frame  overall  depth 

MF      -  Frame  Factor 

NSR    -  Number  of  increments  of  slenderness  ratio 

P      -  Poisson's  Ratio 

PBRY    -  Collapse  pressure  due  to  Bryant's  Formula 

PC      -  Required  collapse  pressure 

PLUNCH  -  Collapse  pressure  due  to  Lunchick 

PTOK   -  Collapse  pressure  due  to  Tokugawa       , 

PWIND   -  Collapse  pressure  due  to  Windenburg 

PY88    -  Collapse  pressure  due  to  #88 

Rl      -  Frame  width/Frame  depth 

R2      -  Web  thickness/Frame  Depth 

R3      -  Flange  thickness/Web  thickness 

SF92A   -  Efficiency  due  to  #92A 

SFWENK  -  Efficiency  due  to  Wenk ' s  Formula 

SIGFOT  -  ap/aT  for  Lunchick  Formula 

SIGBOT  -  obt/^T  for  Lunchick  Formula 

SI92    -  \}j   due  to  #92 

SI92A   -  t[!  due  to  #92A 

SIWENK  -  ^  due  to  Wenk 

SR     -  Slenderness  Ratio 

T      -  Thickness  of  shell 

T2      -  Thickness  of  flange 

W      -  Width  of  web 

WT     -  Weight  of  vessel  per  foot  length 

VOL    -  Volume  displaced  per  foot  length 

YENTRD  -  Subprogram  to  find  neutral  axis 

YMOM   -  Subprogram  to  find  I 
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APPENDIX  G 
COMPUTED  DATA 
and 
EXPERIMENTAL  DATA  FROM  THE  LITERATURE 
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TAB3   >  OF  MAXIMUM  VALUES  FROM  FORMULA  4!9  2A  AND  MAXIMUM 
PR1        S1   ESS  CRITERION 


10  0h/Dm     Tj   max/1. 2   X 


Y' 


Frame  Factor  4 


Frame  Factor  5 


Frame  Factor  6 


251 
313 
,376 
439 
503 
566 
629 
692 
756 
819 
883 
946 
010 

251 
313 
376 
439 
503 
566 
629 
692 
756 

251 
313 
376 
439 
503 
566 
629 
692 
756 
819 
883 
946 
010 


.9086 

.80 

.9490 

.72 

.9806 

.66 

1.0052 

.61 

1.0253 

.58 

1.0426 

.54 

1.0570 

.51 

1.0685 

.49 

1.0800 

.46 

1.0901 

.45 

1.0987 

.44 

1.1074 

.42 

1.1146 

.40 

.9317 

.80 

.9706 

.72 

.9994 

.66 

1.0224 

.61 

1.0411 

.57 

1.0570 

.54 

1.0699 

.51 

1.0814 

.49 

1.0916 

.47 

.9504 

.82 

.9864 

.72 

1.0138 

.65 

1.0354 

.61 

1.0541 

.57 

1.0685 

.54 

1.0814 

.51 

1.0916 

.49 

1.1016 

.47 

1.1102 

.45 

1.1175 

.43 

1.1246 

.42 

1.1304 

.40 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1, 

1 

1, 

1, 

1, 

1, 

1. 

1, 

1 

1, 

1, 

1, 

1, 

1 

1 


785 
664 
574 
509 
447 
410 
375 
341 
321 
291 
266 
252 
240 

716 

602 
519 
458 
409 
367 
334 
304 
278 


642 
554 
485 
419 
373 
334 
305 
1.276 
1.253 
1.234 
1.219 
1.201 
1.191 
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TA1?T  [    MAXIMUM    VALUES    PROM    WENKS    FORMULA 

henc  y-voK  eeld"  crite~rion-"midp"lane~,  midbay 


100h/Dm  ^  max/1 .2       \ 


Y' 


Frame 


Frame 


Frame 


Factor  4 

.376 

1.1333 

.82 

1.646 

.439 

1.164 

.76 

1.611 

.503 

1.198 

.72 

1.562 

.566 

1.214 

.67 

1.533 

.629 

1.237 

.63 

1.506 

.692 

1.248 

.60 

1.481 

.756 

1.265 

.58 

1.454 

.819 

1.270 

.55 

1.438 

.883 

1.280 

.52 

1.425 

.946 

1.292 

.50 

1.409 

1.010 

1.303 

.49 

1.391 

Factor  5 

.376 

1.1592 

.84 

1.620 

.439 

1.1880 

.76 

1.571 

.503 

1.2139 

.70 

1.538 

.566 

1.2341 

.67 

1.497 

.629 

1.2500 

.63 

1.472 

.692 

1.2644 

.60 

1.448 

.756 

1.2759 

.57 

1.428 

Factor  6 

.376 

1.1794 

.84 

1.573 

.439 

1.2068 

.76 

1.540 

.503 

1.2298 

.71 

1.502 

.566 

1.2471 

.66 

1.476 

.629 

1.2609 

.62 

1.451 

.692 

1.2744 

.60 

1.423 

.756 

1.2859 

.57 

1.405 

.819 

1.2960 

.54 

1.390 

.883 

1.3032 

.52 

1.374 

.946 

1.3118 

.51 

1.356 

1.010 

1.3176 

.49 

1.345 

-99- 

FRAME  PARAMETERS  AT  MAXIMUM  EFFICIENCY  #92A 

100h/Dm  Il/IF         lOOb/D 


Frame 


Frame 


actor  4 

.251 

■  .80 

4.758 

.2195 

.313 

.72 

5.335 

.2350 

.376 

.66 

5.860 

.2485 

.439 

.61 

6.340 

.2640 

.503 

.58 

6.817 

.2720 

.566 

.54 

7.238 

.2780 

.629 

.51 

7.655 

.2850 

.692 

.49 

8.072 

.2930 

.756 

.46 

8.464 

.2960 

.819 

.45 

8.871 

.3045 

.883 

.44 

9.267 

.3120 

.946 

.42 

9.662 

.3145 

1.010 

.40 

10.065 

.3150 

Factor  5 

.251 

.80 

5.064 

.2065 

.313 

.72 

5.672 

.2215 

.376 

.66 

6.224 

.2340 

.439 

.61 

6.735 

.2420 

.503 

.57 

7.212 

.2525 

.566 

.54 

7.679 

.2610 

.629 

.51 

8.124 

.2675 

.692 

.49 

8.567 

.2745 

.756 

.47 

9.001 

.2810 

Factor  6 

.251 

.82 

5.364 

.1990 

.313 

.72 

5.958 

.2110 

.376 

.66 

6.533 

.2220 

.439 

.61 

7.062 

.2315 

.503 

.57 

7.566 

.2400 

.566 

.54 

8.056 

.2475 

.629 

.51 

8.527 

.2535 

.692 

.49 

8.995 

.2605 

.756 

.47 

9.451 

.2660 

.819 

.45 

9.919 

.2700 

.883 

.43 

10.391 

.2725 

.946 

.42 

10.853 

.2785 

1.010 

.40 

11.354 

.2790 

-100- 

FRAME  PARAMETERS  AT  MAXIMUM  EFFICIENCY   #WENK 


lOOb/D 


Frame  Factor  4     .376      .82       6.192        .2800 

.2925 
.3055 
.3140 
.3210 
.3290 
.3375 
.3420 
.3440 
.3495 
.3555 

.3        .84       6.5  .2670 

.2755 
.2830 
.2950 
.3010 
.3090 
.3140 

.2540 
.2620 
.2715 
.2780 
.2835 
.2930 
.2980 
.3010 
.3060 
.3130 
.3160 


10  0h/Dm 

li/Ip 

.376 

.82 

6.192 

.439 

.76 

6.647 

.503 

.72 

7.085 

.566 

.67 

7.468 

.629 

.63 

7.843 

.692 

.60 

8.213 

.756 

.58 

8.580 

.819 

.55 

8.930 

.883 

.52 

9.280 

.946 

.50- 

9.634 

L.010 

.49 

9.986 

.376 

.84 

6.568 

.439 

.76 

7.008 

.503 

.70 

7.433 

.566 

.67 

7.878 

.629 

.63 

8.264 

.692 

.60 

8.657 

.756 

.57 

9.041 

.376 

.84 

6.853 

.439 

.76 

7.311 

.503 

.71 

7.768 

.566 

.66 

8.195 

.629 

.62 

8.616 

.692 

.60 

9.037 

.756 

.57 

9.446 

.819 

.54 

9.859 

.883 

.52 

10.273 

.946 

.51 

10.678 

..010 

.49 

11.106 

-101- 
EXPERIMENTAL    DATA    FROM    REFERENCE     9 


Model   Number        lOOh/D      tfy/1000  \         V/92A  Vl 


exp 


1 

.88 

46.8 

.41 

1.43 

1.68 

2 

.784 

47.5 

.42 

1.42 

1.74 

BR7m 

.784 

59.2 

.46 

1.44 

1.62 

BR7 

.599 

44.8 

.59 

1.11 

1.3 

5 

.814 

47.1 

-       .61 

1.01 

1.24 

6 

.474 

49.6 

.67 

1.09 

1.24 

7 

.581 

47.0 

.70 

.997 

1.11 

V7 exp  represents  V/  determined  using  actual  collapse 
pressure  and  V/go*  was  determined  using  collapse  pressure 
calculated  by  formula  #92A. 


-102- 

EXPERIMENTAL  DATA  FROM  REFERENCE  19 

General  Model  Parameters 

Left/Dm     =  .1405 

E  =   30,000,000  lb/in2 

Light  Frames 

MODEL  NUMBER   lOOh/D    tfy/1000  X         Y'exp 


79 

.301 

27 

82 

.303 

28 

76 

.307 

42 

88 

.293 

27 

75 

.294 

39 

81 

.295 

28 

84 

.300 

29 

87 

.300 

29 

80 

.310 

28 

85 

.311 

29 

92 

.288 

32 

91 

.290 

32 

93 

.308 

29 

96 

.298 

28 

94 

.298 

28 

90 

.300 

27 

876 

.92 

855 

.90 

06 

.66 

890 

.85 

06 

.91 

900 

.93 

909 

.81 

905 

.81 

840 

.87 

884 

.83 

00 

.85 

975 

.84 

886 

.88 

896 

.85 

896 

.84 

88 

.82 

In  the  computation  of  Xvalues>  n/D  was  used  instead 

of  h/Dm  resulting  in  a  small  error  which  was  not  critical 

in  the  graphical  methods  used. 

V^xd  rePresents  tne  ratio  of  actual  collapse  pressure  to 
the  hoop  stress  collapse  pressure. 


-103- 

ALUM INAUT 

In  the  following  section,  X,  ty ,    and  n  have 
been  found  for  the  pressure  hull  in  Aluminaut.   The 
dimensions  of  that  craft  were  taken  from  References 
(5)  and  (15)  and  are  for  the  cylindrical  hull  alone, 


Operating  Depth  =   15,000  ft.       LOA 
Safety  Factor    =1.4  Le 

av  =   60,000  psi      D 

00  psi   Dm 
3       h 


E 

=      2 

9,000, 

P 

X 

1   lb/ii 

1. 

4' 

(L/Dn)2     1 

°Y 

\ 

(VDm)3      \ 

E 

X 

= 

.324 

2. 

* 

= 

Pc 

0h 
2D°Y 

1 

] 

33.3  ft, 

6.6  ft, 

8.0  ft, 

7.5  ft, 

.5  ft, 


p   =  1.4  X  15,000  X  .446 


pc  =  9400  psi 


ty      =    1.25 


3.  n  = 


PcV 


W  10  - 

r)      =        2.65  in. 


From  the  above  computations,  it  is  apparent  that  the 
designers  of  Aluminaut  produced  a  vessel  whose  parameters 
put  it  in  the  low  lambda  design  region,   However,  it  has 
been  noted  (15)  that  a  more  efficient  structure  is  possible, 
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